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ABSTRACT 


In this dissertation five simulation technigues ((1) 
Sampling without replacement from the observed distribution 
of impact energy densities; (2) sampling with replacement 
from the observed distribution; (3) sampling from aA 
idealized two-parameter gamma distribution; (4) stratified 
sampling from the observed distribution of impact energy 
densities; and (5) spatial sampling from a contoured pattern 
of the observed impact energy density field) are presented 
and applied using the Monte Carlo method in an attempt to 
simulate the sort of sampling a hailpad network would 
produce, Data from the 29 July 1974 hailswath is used to 
illustrate the approach for each of these five techniques. 
An estimate of the errors involved in the mean impact energy 
densities for the hailpad observations is also incluied 
using the first four techniques. Data from five hailswaths 
that occurred in central Alberta during the three years 1974 
to 1976 are used in the analysis. 

The collection and calculation of the impact energy 
densities from the hailpad data are also described. Graphs 
of the probability density function of the sample means, the 
probability distribution function of the sample means, and 
the standard deviation of the sample means as a function of 
the sample size are given for each simulation technique 


outlined above. 
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The results show that the probability density function 
of the sample means approximates the Gaussian form {in most 
cases) for sample sizes greater than about 15. The results 
further show that the variation of the standard deviation as 


-Y; 
* for 


a function of the sample size is proportional to k 
nhon-stratified sampling, stratified sampling, and sampling 
from the fitted two-parameter gamma distribution. These 
results show good agreement with the predictions of the 
central-limit theorem. 

Toa test how adequately the two-parameter gamma 
distribution fitted the original data, the Chi-square test 
was used. The results for the five hailswaths show that this 
distribution gives a reasonable fit for the impact energy 
densities. 

Which technigue best simulates the sort of sampling 
which hailpads do in a hailswath, depends on the nature of 
the hailswath itself. Sampling with replacement (non- 
stratified), is probably best applicable to hailswaths where 
the impact energy density distribution in any small sub- 
section of the hailswath is essentially the same as that in 
the entire hailswath. Sampling from a fitted two-parameter 
gamma distribution is probably best applicable to severe 
hailswaths where a large range of impact energy densities is 
observed. Stratified sampling can be ae as a Simulation 
technique for hailswaths where the impact energy densities 
are continuous in space. In deciding which technique best 


Simulates the hailswath it is necessary to consider the 
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auto-correlation function for the impact energy densities. 
Using the Monte Carlo method, the results indicate that 
in order to commit a 10 percent error between the average of 
the sample means and the the mean of the sample means. 90 
hailpads are required in a 1000 km? area (or 1 hailpad per 
11.1 km2) for non-stratified sampling and 200 hailpads per 
1000 km2 area (1 hailpad per 5.0 km?) are needed for the 
fitted two-parameter gamma distribution. An improvement in 
the sampling density is observed for stratified sampling, 
where for a 10 percent error between the estimates, 10 
hailpads are needed in a 1000 km? area (1 hailpad per 100 
km2). The major conclusion of this dissertation is that the 
present dense network (with 1 observation per 5.7 km2) of 
hailpads is sufficient in density to produce an error of 
less than 10 percent in the impact energy density. The 
precipitation network (of t observation per 70 km?) is, on 
the other hand, insufficent in density at this error level 
(for sampling from the fitted two-parameter gamma 
distribution and for non-stratified sampling). Assuming a 
100 percent response in the hailpad observations for the 
area affected by the swath, the minimum error in the impact 


energy density with the present precipitation network would 


be at least 18 percent. 
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CHAPTER 1 


INTRODUCTION 


1.1 Background 

Hail has been reported in central Alberta during the 
last 20 years on an average of 60.2 days 1 each year between 
15 May and 15 September (Wojtiw,1977). Although hail damage 
on many of these days is of no consequence, claims for hail 
damage in the same period are received by the Alberta Hail 
and Crop Insurance Corporation on an average of 50 days each 
season throughout Alberta, with the majority of these hail 
days occurring in central Alberta (Summers and Wojtiw, 1971). 
The total annual loss to agriculture from hail is now 
estimated to average between 50 and 60 million dollars each 
year (Wojtiw,1975b). Concern over such heavy agricultural 
losses and an enthusiastic meteorological research interest 
in severe storms led to the creation in 1974 of a hail 
suppression and hail research program called the Alberta 
Hail Project (AHP). The program, supported by the Government 
of Alberta, is under the auspices of the Alberta Weather 


Modification Board (AWMB) with participation by the Alberta 


1 The standard deviation is 9.1 days, while the median value 
is 58 days. 
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Research Council({ARC), the Atmospheric Environment Service 
of Environment Canada {AES), and the Universities of Alberta 
and McGill. This program was created from the Alberta Hail 
Studies(1956-1973) program(ALHAS), and the research programs 
in existence in 1973 were continued in the AHP program. The 
ALHAS program was directed by the Alberta Research Council 
with financial and scientific contributions supplied by the 
Atmospheric Environment Service of Environment Canada, the 
National Research Council of Canada and the Stormy Weather 
Group of McGill University. 

The AWMB was formed in March 1973 by the Legislative 
Assembly of Alberta, in response to the report of the 
Special Legislative Committee on Crop Insurance and Weather 
Modification (Stromberg,1972). The Board was established 
under the Alberta Department of Agriculture with a five year 
mandate to administer hail suppression and hail research 
programs. A truncated circle containing 47,900 km@ (Figure 
1.1), centered on the meteorological radar site at the Red 
Deer Industrial Airport was designated as the project 
boundary. The circular area was divided into a northern and 
a southern part with a buffer zone between. Each zone has 
different cloud seeding criteria and rules (Deibert and 
Renick,1975; Deibert,1976; Deibert,1977). 

The Alberta Hail Project collects various data (radar 
and rawinsonde measurements, hailpads, hailcards, among 
others) to further increase knowledge of hailstorms and to 


evaluate the success of the hail suppression technique. Of 
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the ground measurements collected, hailpad data (see Chapter 
2 for details) are considered to be the prime source of 
objective hailfall data, and they are to be used in the 
evaluation of the hail suppression hypothesis. The hailpad 
data to be used in the evaluation come from a network with 
an average spacing of one hailpad per 70 km2. Before the 
data can be used in any evaluation scheme, it is desirable 
to know if there are sufficient hailpads in the network to 
measure adequately the hailfall and what degree of error is 
associated with the hailpad values. The answer depends on 
what is meant by "adequately". If, for instance, it is 
desirable to compare areal energy density patterns on the 
ground (from the hailpads) with radar echoes aloft, the 
spatial resolution of the hailpad network should be as good 
as that of the radar. The radar # echoes aloft are recorded 
by a 10 cm radar with a resolution of 1.6 km@ at a range of 
50 km, while the areal energy density of the hailswath is 
estimated using data from a network with an average spacing 
of one hailpad per 70 km2. Since the average spacing of the 
two sets of measurements differs by at least an order of 
magnitude, they are not readily comparable. If, on the other 
hand, it is desirable to obtain an accurate estimate of the 
mean areal energy density of the haiiswath, the required 


2 Characteristics of the S-band radar: frequency=2.880 GHz; 
pulse repetition frequency=480 s~’ ; pulse duration=1.754,s; 
beam width=1.15° (all planes); antenna rotation rate=8 RPM; 
and elevation program-spiral scan of 1 per revolution to 8° 
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Figure 1.1 Map showing area of the Alberta Hail Project 


network density can be specified in terms of the inaccuracy 


that one is prepared to tolerate. 
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1.2 Purpose of Study 

It is this latter aspect, namely, the desire to obtain 
an accurate estimate of the mean areal energy density of the 
hailswath and the density required, to do this which is 
addressed in this dissertation. Here five different 
Simulation techniques are presented and applied using the 
Monte Carlo method (discussed in Chapter 3). The purpose is 
to estimate the errors involved in determining the mean 
areal energy density of the hailswath from hailpad records 
With a given areal density. The collection and reduction of 
the data used in this dissertation are described in Chapter 
2 while in Chapter 4 the five techniques are examined using 
data from the 29 July 1974 hailswath as an example of the 
approach. The results from the analysis of five hailstorms 
(1974-1975) are given in Chapter 5, with the conclusions in 
Chapter 6. 

This work is a first effort, as far as can be 
determined, in using the Monte Carlo method to examine the 
sampling density problem for hailpad data. The technique is 
not new; Similar ideas have been used elsewhere, for example 
in theoretical work on radar signal statistics (Marshall and 
Hitschfeld,1953; Smith, 1964); comparison of rainfall 
measurements on the ground with radar signals 
(Collier,et.al.,1975), and with hail crop damage data 
(Morgan and Towery,1974). To assist in obtaining a finer 
scale picture of the areal energy density distribution than 


is permitted by the large scale network, data from a dense 
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network(see Figure 2.2) have been collected during the years 
1974 and 1975. The average spacing in this network is one 
hailpad per 5.7 km2. 

The problem of assessing the accuracy of hailfall 
networks has received some attention in the United States 
and in Alberta, Canada. A brief summary of this work, 
including a discussion of the hailfall phenomenon, is given 
in Section 1.3. Haiifalls are believed by various authors 
(Changnon, 1968b; Morgan and Towery,1974), to vary more in 
space and time than do rainfalls. Due to this variability, 
assessing the accuracy of a network is a more difficult 
problem for hailfall than for rainfall; thus, the limited 
efforts in this area of research to date. Considerable 
effort has been devoted to the problem of assessing the 
accuracy of measurements obtained using raingauge networks. 
A review of the literature on this subject is included in 
Section 1.4 in order to give an idea of the type of 


behaviour which may be expected for hailfall. 


1.3 Precipitation Patterns - Hailfall 
Extensive research programs into the climatology of 


hail in Alberta 3 , Colorado * , and Illinois S$ have been 


—oo 


3 (Douglas and Hitschfeld,1959; Longley and Thompson,1965; 
Paul, 1967; Summers and Wojtiw, 1971; Wojtiw,1975) 

* (Schieusener,1963; Renne and Sinclair,1968; Beckwith,1960; 
Swinbank,1971; Schickedanz and Changnon, 1970, 1971) 

Ss (Huff and Changnon, 1959; Huff,1960,1961; 
Changnon, 1962, 1967, 196 9a, 1971) 
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pursued continuously since the late 1950's, and they have 
produced much of the total climatic information about hail 
available in the world. These hail studies have yielded 
extensive, largely unsummarized, information relevant to the 
time and space variations of hail. 

One of the first programs to aim at eventual hail 
suppression efforts was begun during the late 1950's in 
Colorado (Schleusener and Jennings, 1960). With it, the 
hailpad (see Chapter 2) as a hail sensor was introduced to 
obtain quantitative measurements of hail size distribution 
and number concentration. In the 1960's more hail networks, 
also composed of hailpad sensors, were instituted in western 
South Dakota (Koscielski,1967), western North Dakota 
(Schleusener, et al.,1965), and Nebraska (Musil and 
Dennis, 1968). More densely spaced networks, largely 
utilizing rain/hail separators and raingauges, were 
installed in central Illinois in 1967 (Changnon,1969a) and 
in north-eastern Colorado in 1969 (Swinbank,1971) to provide 
the basic data for the evaluation of suppression 
experiments. 

In the analysis of 17 years of volunteer reports of 
hailfall, Summers and Wojtiw (1971) found that hailswaths 
which caused damage to crops in Alberta were most frequently 
from 3 to 8 km wide, with some severe storms producing 
swaths with crop damage 30 km wide. A hailswath has been 
defined as a precipitation pattern observed on the surface, 


resulting from convective activity from which hail was 
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reported (Wojtiw,1975a). Wojtiw used six hail reports from 
the same storm system as a minimum criterion to define a 
hailswath. Areas of 3 to 10 km2 where rain without hail can 
exist, have been observed in the hailswaths of some severe 
Storms (Pell,1969). Within many hailswaths, Changnon, et.al. 
(1967), using very dense instrumented networks in Illinois, 
identified a hail entity called a "hailstreak", and defined 
this as "the surface outline of a single volume of hail 
produced in a storm". Figure 1.2 (Changnon,1968d) shows an 
example of hailstreaks (small enclosed areas) in two major 
hailswaths (the two large strips extending across the 
figure) indicating that the two large storms which traversed 
this 57,000 km2 area of Illinois, produced semi-continuously 
several separate volumes of hail. Also shown in the figure 
are many other hailstreaks which are associated with other 
systems which crossed the region in the 10-hour period. The 
region between the hailstreaks and the envelope of the 
hailfall pattern are undefined by the author, but probably 
consist of areas of hail or rain used in defining the 
envelope of the hailfall pattern. The hailstreaks shown in 
Figure 1.2 were drawn using data from a dense network of 
hailpads (one hailpad per 7 km2), supplemented with 
insurance loss and field survey data. From observations of 
434 hailstreaks in central Illinois and around he St.Louis 
area in 1967-1968, Changnon (1970) computed that "their 
median area was 10.5 km2 (ranging from 2.3 to 2040.9 km?), 


with an average width of 2.8 km and an average length of 
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Figure 1.2 Diagram showing hailstreaks and hailswaths in a 
10-hour period on 21 April 1967 in a region of Illinois. 
Hailstreaks are the small enclosed patches; hailswaths are 
large strips that traverse the map. (after Changnon, 1968d) 


15.3 km". Further time-space analysis of hailstreaks led to 
the definition and study of "hail areas" 

(Changnon, 1968c),that is, the surface outline of the 
hailstreak at any given instant. Admirat (1972) has applied 
the term "hail core" to the concentrated high-loss (severe 


damage) segment of the hailstreak. 


An even finer areal study of hailfall patterns (Morgan 
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and Towery,1974) was undertaken using very dense networks 
(hailpads spaced 10 to 100 m apart) over 0.3 to 3 km2 areas. 
Networks of this type in Colorado and Tilinois sensed only 
Segments of a hailstreak which crossed the networks, but 
they provided extremely detailed information on the 
structure of streaks. Figure 1.3 shows the hailstone number 
density pattern for a hailstreak that crossed a 1.6 by 1 km 
network with 120 hailpads, and it depicts the considerable 
variability in the frequency of stone occurrence (values 
from 1076 to 3767 m*) ©. 

In Alberta, limited efforts in assessing the sampling 
density have been conducted (Strong,1974). This dissertation 
is a continuation of that work, and is a serious first 


effort to determine the necessary sampling density for 


hailpad data. 


1.4 Precipitation Patterns - Rainfall 

In this section, a review is given of the literature on 
network measurements of rainfall, as the results may provide 
an idea of what can be expected for hailfall. For rainfall, 
a number of studies have addressed themselves to problems of 
the relationship between the network density (nuaber of 
observations per unit area) and the representativeness of 


the mean value of the rainfall from a particular event fora 


6 Values on Figure 1.3 are given in? und bsiipe rm: it 4, ache 
conversion is 100 ft'*=1076 mm“. 
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Figure 1.3 Diagram showing the number of measurable stones 
from 3 April 1974 storm «after Morgan and Towery,1974) 


given area. The required density cf a network for rainfall 
is determined mainly by the variability of the precipitation 
events and the measurement objectives. Generally long- 
duration and ee eeneea precipitation is much less variable 
than short-duration and small-area precipitation; therefore, 
the density requirements for the latter case are much 


higher. For example, it was found for the 2927 km@ (208 mi#) 
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Chickasha watershed that at the 5 percent significance 
level, 10 gauges were as good as 158 in assessing mean-daily 
areal rainfall. However, for other objectives such as 
Sampling peak storm-rainfall, the reduced network of 10 
gauges was of questionable value (Nicks, 1965). 

Neff (1965), using network density of one recording 
gauge for each of the 241 km2 of area in the Reynolds Creek 
basin, found that as the gauge density is decreased, the 
accuracy of estimating the areal mean is likewise decreased. 
The estimates Neff obtained had the 95 percent confidence 


limits given in Table 1.1. These results show that the 95 


Table 1.1 Estimates of the Error of {or Uncertainty of) 
Areal Mean Rainfall for a Given Number of Gauges 
No. of Gauges 95% Confidence Limits 
90 + 2%90 ¢D 
30 + /3.6/ cn 
10 + 6.9 co 
4 $13.2 0n 


percent confidence limits in the mean areal rainfall are 
inversely proportional to N’’, where N is the number of 
raingauges. 

Numerous empirical studies of dense networks have been 
made in order to determine the number of raingauges 
necessary to estimate areal rainfall amounts with a specific 


accuracy. From one such study of 15 networks in the United 
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States, Hershfield (1965) adopted an inter-gauge correlation 
coefficient of 0.9 as the criterion for raingauge spacing. 


His results are shown in Figure 1.4. 
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Figure 1.4 Diagram showing results of estimating required 
distances between gauges as a function of the 2-year-24 hour 
and 2-year-1-hour rainfall (Note 0.9 inter-gauge correlation 

coefficient) (after Hershfield, 1965) 
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Hutt 6(197-0),,. using data from two dense networks of 
raingauges in Illinois, obtained estimates of sampling 
errors in the measurement of areal mean precipitation on 
areas of 130 to 1295 km2 for daily, monthly and seasonal 
precipitation. The storm sampling error was found to 
increase with increasing areal mean precipitation and to 
decrease with increasing sampling density (gauges per unit 
area) and increasing storm duration (Figure 1.5). In this 


figure the sampling errors were fitted using Equation 1.1. 


log E=-1.5069+0.65P+0.82G-0. 22T-0.45A (1.1) 


where A is the area (mi?), E is the average sampling error 
in inches, P the areal mean precipitation in inches, G the 
gauge density in mi@ per gauge, and T the storm duration in 
hours. Equation 1.1 was computed by combining all network 
data from the East Central Illinois Network (49 recording 
raingauges in 1036 km@ (400 mi#)) in the May to September 
period. Relatively large differences were frequently found 
in the sampling errors among storms of apparently similar 
characteristics (i.e. frontal storms, air mass storms, and 
low pressure center passages). Huff found that air mass 
storms require the greatest sampling density among synoptic 
storm types to maintain a given error level. In the warm 
season, rainshowers and thunderstorms were found to require 
nearly twice as many gauges as steady rain for a given 


measurement reliability. Also, the sampling density 
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Average Sampling Error, Inches | 


AOD se OO I OMG) a fee: 


_ Storm Mean Rainfall, Inches 


Figure 1.5 Graph showing sampling error relations on 1036 
km2 qafter Huff,1970) 


requirements for the May to September period were two to 
three times those needed in the October to April period.. 

A considerable difference in the magnitude of storm 
sampling errors during the May to September period was found 


between consecutive 5-year periods on the same network. 
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Analyses indicated that this difference resulted partially 
from a much higher number of air mass storms in one of the 
periods. Therefore, caution should be used when interpreting 
sampling error relations based upon observational periods of 
five years or less. Huff also concluded that sampling 
requirements for the measurement of monthly and seasonal 
rainfall are less stringent than those for individual storm 
precipitation. He found that for seasonal rainfall and a 
relatively sparse density of 518 km? per gauge, the average 
sampling error was less than 5 percent for the entire range 
of seasonal rainfalls. By way of contrast, for a 6 hour 
rainfall with a mean of 1.27 cm, the average error with the 
same network density was 14 percent. Huff mentions that his 
results apply to the establishment of raingauge networks in 
the Midwest United States and in other areas which have a 
Similar precipitation climate. 

In more recent studies, Eddy (1976), using rainfall 
observations from the Illinois State Water Survey at Goose 
Creek and the East Central Illinois networks, devised a 
model using optimal interpolation objective ee in 
order to estimate the rainfall, the gauge density and the 
accumulation time. In this model, the a priori rainfall was 
modelled by fitting a gamma distribution function to the 
observed rainfall amounts from individual storms and to the 


composite data from a combination of storms. The Type 1 and 


yilsisisq botivas2 ssaszettib 


quod 38.5 tot sidan daira 40> Bie a 2 


ods dtiw torts spaieve aas a3. s.r 


ain t5d3 enotinen awh  tng075 


e200 +6 yovaue 49350 va atoatir wae wa ‘seotsevsondo, 
s bealveb ,edzewsen etoakitt 1 ieasas> sess os bas 499729 7 a 
me sieabaihsst ovivootdo a pf ai ee gies. Li Lesisqo Rate fefom uf 
out bas ysienesbh semen oe Mila cd osenises oF 296% 
oboa aids ar “sont noitslumypot ee) i 
5 ices 6 patssid vd ba Lisbon 2 a 


5 | azavoas listaisa bovxead 


asW iteiates rioiag * edt a 


bins u ‘ edz oh items? soitedias 
ie B iW 2 ‘ae * 


ots OF bas anio%2 sevnaviba 


“ mo3s ial: f 
va fh 5 , f na 


17 


Type 2 errors 7 are graphically depicted in Figure 1.6. He 


RAINFALL FREQUENCY 


Figure 1.6 Distribution of two objective analysis sets about 
their respective true values with a priori true rainfall 
distribution superimposed. The Type 1 («) and Type 2 2 ) 
errors are indicated. {after Eddy,1976) 


concludes that for a dense network (816 km2 area) the 
optimal values ® give an average density of 16.3 km2@ per 
gauge with an accumulation time of about 10 min while fora 
sparse network (25,900 km2 area) one gauge per 104 km? will 


give the desired accuracy % if at least one hour of rain is 


7 “Given an acceptable objective analysis R e one can say 
that it represents a random deviation from population R, if 
R > Re and from Rk, if 8.2 °R- with the following expectation 
of misclassification: The fraction of the time R, would be 
chosen (using Rap R.), when R, was the truth, would be 

ec (Type 1 error). The fraction of the time when R, would be 
chosen (using R < R,), when R,z was the truth, would be 

Bo (Type (ZeeLrrorjy.” 
8 Minimum rainfall of interest Rm=0.02 inches; Type 1 error 
(«< =0.01), Type 2 error ( ps =0.10) 

9 Minimum rainfall of interest Rm=0.1 inches; Type 1 error 
(eats 1h) Ty pegeeerror a(¥ 2 =0.'15) 
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collected by any gauge during a 24 hour period. The analyzed 
rain came from unstable summer showers and was observed 
using tipping bucket raingauges. Consequently, Eddy's 
humerical results apply to shower-type rainfall and should 
not be transferred to other locations without a careful 
consideration of whether the showers have the same 
statistical structure as that found in Eddy's data sample. 
In another study, Collier, et.al. (1975) compared 
calibrated radar measurements of areal rainfall with 
rainfall measurements obtained solely with raingauge 
networks of various densities. In that study, the authors 
compared the accuracy of the two methods of measurement in 
hilly terrain over a period of one hour, using results 
obtained with a one degree beamwidth and a 5.6 cm wavelength 
radar. The comparisons were made over areas which varied in 
size from 20 to 100 km?@. The raingauge network was 
distributed over 1000 km2 and consisted of 62 battery 
operated Plessey MM37 tipping bucket raingauges, 
incorporating a magnetic tape event recorder. Figure 1.7 
shows the results of this study; that is, the accuracy of 
the investigations achieved by a calibrated radar 19 and by 
raingauge networks of various densities. The curved lines 
(solid indicate the calibrated radar, and hyphenated signify 


the raingauges) show the mean percentage differences from 


10 The radar calibration is performed with one or more rain 
gauges. 
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Meon/Difference/between estimates 


1 2 5 10 20 50 


Number of raingouges over 1000 km? 


Figure 1.7 Graph showing difference between estimates of 
hourly sub-catchment rainfall related to raingauge density 
(after Collier,et.al., 1975) 


the optimum estimates 4! of hourly sub-catchment !2 rainfall 


as measured by various numbers of raingauges distributed 


over the 1000 km? area, and the calibrated radar system. The 


results from the raingauge networks show a large difference 


in accuracy between days with showers and those with 


11 The observed radar echo intensity pattern was used to 
interpolate between the gauges so as to obtain a reliable 
estimate of rainfall features on a scale smaller than the 


gauge spacing. The resultant distribution of precipitation 
while the values in the field 


are the optimum estimates. 
12 one of 16 unequal segments of the 1000 km? area. 
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widespread rain. The accuracy which could be achieved by a 
radar calibrated by raingauges and incorporating a model 
(Harrold and Kitchingman, 1975) to reduce the melting layer 
error is shown by the upper solid curve in Figure 1.7. 

The solid lines on the diagram are the mean percentage 
differences between optimum estimates of hourly rainfall and 
radar estimates calibrated using various numbers of 
calibration sites. The higher accuracy line (lower solid 
line on diagram) indicates occasions when the radar beam is 
entirely within rain. The lower accuracy line (upper solid 
line on diagram) indicates occasions when a mixture of rain 
and snow (sleet) is within the bean. 

The hyphenated lines on the diagram represent the mean 
percentage differences between optimum estimates of hourly 
rainfall and estimates using various numbers of raingauges 
without a radar. Upper lines show showery conditions, lower 
lines indicate widespread rain. The diagram (Figure 1.7) can 
be used for example to derive the number of raingauges 
required to measure rainfall with the same accuracy as a 
radar system calibrated with a specified number of raingauge 
sites. For example in typical showers a radar system 
calibrated using two raingauge sites gives the same accuracy 
as a raingauge network of about 50 raingauges over 1000 km@. 

In this chapter several research projects dealing with 
precipitation patterns have been summarized. The efforts 
made by various researchers to assess the relationship 


between hailfall and rainfall variability, and network 
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density have also been presented. The limited effort in 
resolving this problem for hailfall has been twofold: 
unavailability of quantitative data and the variability in 
Space and time of hailfall. Much more work has been done on 
assessing the required network density for rainfall. Before 
introducing the approach for assessing the hailpad network 
density, a discussion of the collection and reduction of the 
hailpad data will be given, so that the reader may become 


familiar with the observations used. 
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CHAPTER 2 
DATA COLLECTION AND REDUCTION 


2e1 Observations 

The primary source of objective hailfall data is the 
hailpad-a styrofoam pad of dimensions 30.5 cm x 30.5 cm x 
2-5 cm. Hailpads were set out in two networks operated by 
the Alberta Hail Project during the summer months in central 
Alberta from 1974 to 1976. 

The two networks and their locations are shown in 
Figure 2.1 and Figure 2.2. These are: 

(1) a precipitation network of 500 stations over an area of 

34,467 km@ 2 
(2) a dense network of 162 stations over an area of 932 km2 

(embedded within the large precipitation network in the 

Rimbey district of Alberta). 

By analysis of the pads, the hail size distribution, 
and the areal concentration of the hailstones for the 
duration of a hailstorm can be obtained. However, the pads 
do not give time-resolved information about the hailfall, 
which may be important for relating cloud seeding events to 


hailfall on the ground. 


All stations (volunteer observers) in the two networks 


1 No stations are located on the western edge of the project 
area, hence the figure reported here is smaller than that 


given in Chapter 1. 
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Figure 2.1 Map showing networks in the Alberta Hail Project 
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Figure 2.2 Diagram showing location of hailpad stations in 


the dense network (Stations are located in the quarter 
section where the numbers are printed) 
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were supplied with four styrofoam pads, a pad holder #2 
mounted on a supporting stand 3 (see Figure 2.3), a 
raingauge * and appropriate forms to report the occurrence 


of hailfall and the total amount of daily precipitation. 


Figure 2.3 A picture of a hailpad stand 


a 


2 A 30.5 cm x 30.5 cm x 2.5 cm plywood board with two narrow 
metal slits on two of the sides of the board to hold the 
hailpad in place. ; 

3 A 5.1 cm x 10.2 cm x 1.8 m two-by-four, or fence pole if 


available. 
*# A plastic wedge type manufactured by Edwards Manufacturing 


Co., Minnesota, U.S.A. called Tru-Check. 
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Each station operator was asked to change the pad 
immediately after the hailstorm or once per week, or if the 
pad became damaged in any way. Field tests in 1975 showed 
that the pad begins to deteriorate after 10 days in 
sunlight. Beyond ten days, the deterioration is rapid and 
the sizes of hail dents made thereafter are unreliable. 

When they arrive from the manufacturer, the styrofoam 
pads are quite coarse. In order that the roughness is not 
mistaken for smali hail dents, a pad sander (Figure 2.4) was 
developed in 1975 at the Alberta Research Council. The 
dented pads, returned from the field, were then inked with a 
print roller to provide contrast between the uninked hail 
dents and the inked undented background styrofoam. Figure 
2.5 shows an example of an uninked and an inked pad. A 
quantitative measure of the number of hailstones in various 
size categories can be obtained from the hail dents. These 
measurements are used to calculate the impact energy per 
unit area (subsequently referred to as the impact energy 
density) (see Section 2.3). 


2.2 Analysis of the Hailpads 

The analysis of the pads for hail size distribution and 
concentration was carried out by hand. This involves the use 
of a transparent overlay (see Figure 2.6) with imprinted 
Circles of diameters 0.25 cm to 2.54 cm at .25 cm intervals, 
and 2.54 cm to 5.08 cm at .63 cm intervals. The dent was 


estimated to the nearest .13 cm and placed in the 
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Figure 2.4 A picture of a pad sander used to eliminate 
coarseness in the styrofoam 


appropriate size category. For example, if a dent was larger 
than .25 cm but smaller than .38 cm, it would be placed in 
the .25 cm category. This procedure gives rise to an 
underestimate of the sizes, and as a result an underestimate 
of the calculated impact energy density. The magnitude of 
the error in the impact energy density decreases as D*. The 
maximum error in the impact energy density occurs if all the 
size ee are underestimated by one size category. In 
this case the error in the impact energy density is 80.2 
percent for the .25 cm category and decreases proportionally 


as D*Tto anserror eof 17.7 fpercent cinnthe 2354 ,cnceategorys 
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Figure 2.5 A picture of an uninked hailpad (above) and the 
same pad after inking (below) 


For the 29 July 1974 hailswath the maximum error in the 
impact energy density, assuming that all the dent size 
categories are underestimated by one size category, is 36.9 


percent. 


The number of dents of each size category was recorded 
for each hailpad. In order to convert from dent sizes to 
hail sizes, the calibration curve shown in Figure 2.7 was 


used (Lozowski and Strong,1977). 


2.3 Calculation of Impact Energy Density 


The kinetic energy of a hailstone upon impact (impact 


SRP sah eft Ils sahy 5 apie ‘sesh pawn sisindt 


@.0€ af ,yiopets. exta sno yd bssemisestobay ers. aefropesso 
a | | | - > 
thd . Live a, 


3 


‘bsbtoost asy y1opesso ote dose to ‘atosh to | 

ot gosta ta9b aoxd $a9v209 o+ 25520 at. beqhisd dobs s0R ~ 
a6W wo spent ni awode” eNabss noitszdiiso odd yada Lied 

s ileclaaeiai bas Edzwosod) tou ) 


Ll 


bp | r pene esos soa to aottselvaisd €.8 


Ns 


OV EUR LA 


Helm nom eaeAeD 


S) 
ro) 
— 


2,00" 


Figure 2.6 Hailpad overlay 
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Figure 2.7 Graph showing conversion from dent diameter(on 
the hailpad) to hail diameter (after Lozowski and 
Strong, 1977) 


energy), assuming no horizontal wind component, can be 
calculated from the hailpad data using the following 
formulae (Strong,1974). For a spherical hailstone of mass a, 
diameter D, density /, and falling with terminal velocity V, 
through air of density Ze the HEC of motion can be 


solved to yield: 
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ta 
VL= (oes (2.1) 
T 3 fab ‘ 


The following parameters were found to be most appropriate 
for typical Alberta conditions (at a pressure of 890mb and a 
temperature of 22.6 ) (Strong,1974): drag coefficient 

Cy =0.60, hailstone density / =0.89 g ca’, surface air 
density / =1.05x x10 gm cm and g=981.3 cm mis Using these 
values, Equation 2.1 becomes: 
D2 


V,= /3.6 


r (2.2) 


where V; is in meters per second and D in centimeters. _ 
This relation was confirmed by Lozowski and Beattie 
(1975), who used high-speed photography to measure the 
terminal velocity of hailstones up to 1 cm in diameter for 
several 1974 Alberta hailstorms. They found that the curve 
le 


(Se ne a ae 5 was a reasonable fit to their measurements. 


a 
The kinetic energy for a single stone on impact is 


therefore: 


I= ( are) IN (2.3) 
é 7 fa Cd 


Substituting the typical values for the various parameters 


into Equation 2.3 yields: 


-2 n# 
te 4.307 xo D (2.4) 
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where J, is in Joules and D in centimeters. 

This equation gives the kinetic energy of a single 
hailstone. In order to calculate the total kinetic energy 
for the stones on the hailpad, for each of the size 
categories, an impact energy is determined for each category 
and these values are then summed. Equation 2.4 for the 
impact energy density on the hailpad then becomes: 


Pp 

4 

t= 4636x10' 2 ALD; (Tms) (2.5) 
cel 


where D; is the diameter, and A; is the number of hailstones 
on the hailpad in the i th size category (for i=1,P where P 
is the number of size categories). 

For a hailstone falling with velocity V,, ina 
horizontal wind V, and assuming that the horizontal velocity 
component of the hailstone equals the wind speed, the total 
aapaee energy becomes: 

Pror@Sa (V2 tay, 2) (2.6) 
3 
where V, = V-tan@ and m = Tee and where @ is the angle to 


the vertical at which the hailstone strikes the horizontal 


pad. The impact energy then is 


I= Gan D7 ac? (2.7) 
*pa“D 


or upon substituting the typical values given previously 


into Equation 2.7, we obtain: 
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33 
T= 0.4636D*(sec29) (J m*) (2.8) 


Equation 2.8 is the impact energy density for an 
individual hailstone falling with velocity V, ina 
horizontal wind V,. To obtain the total impact energy 
density for the hailpad, the value obtained in Equation 2.8 
is summed over all hailstones and their angles of impact. 
Equation 2.9 is the equivalent counterpart to Equation 2.5 
for a haiistone falling in a wind: 


p 
Es, ie Fo. 
Teng 463610 2 ALD: a0? O (2.9) 


cat 


where’ D7 is\the diameter, A; is the number of -hailstones 
in the ith size category (for i=1,P where P is the number 
of size categories). 

The equations for a hailstone falling in a horizontal 
wind ~ given for completeness, since these have not yet 
been incorporated into the calculation of impact energy 
density from the hailpads. Hence the impact energy density 
values used in the following chapters will be values 
calculated assuming no horizontal wind component. The reason 
for this lack of implementation is that few hailpads have a 
wind measurement associated with them. Wind speed, however, 
is an important quantity. Strong winds an cause 30 percent 
more damage to crops than light to moderate winds (Wojtiw 
and Renick,1973). Morgan (1976) suggests that wind 


contributes to even greater crop damage but the magnitude of 
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the effect is not given. Hence this quantity should not be 
ignored in the evaluation of hail suppression. Using a 
hailcube, five hailpads joined to form a cube (minus the 
bottom hailpad), Morgan has been able to estimate the 
horizontal wind component from the hailpads {Vento and 
Morgan, 1976). These cubical units have not been used in 
Alberta in the networks but should be considered in setting 
up future networks. Hence the (vertical) impact energy 
values calculated by Equation 2.5 are used in the Monte 


Carlo method described in the next chapter. 
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CHAPTER 3 


THE MONTE CARLO METHOD 

3.1 Introduction 

Hailpad measurements are a "random" sample of energy 
densities from a true distribution of energy densities which 
could be determined in principle only if the ground were 
completely covered with hailpads. From the point of view of 
comparing seeded and unseeded hailstorms, a true mean impact 
energy density for the hailswath is desired. Since the true 
mean impact energy density is not available in reality, a 
procedure like the Monte Carlo method is used to obtain an 
estimate of the uncertainty of the mean impact energy for 
the hailswath. This is done by averaging the hailpad 
observations, perhaps weighting them according to the area 
they represent, and estimating how the hailpad average 
differs from the true mean. In this dissertation a certain 
true distribution function is assumed (obtained from the 
hailpad observations) and the Monte Carlo method is used to 
simulate the "random" sampling. Thus for a given true mean 
value a distribution of observed means is obtained. This 
distribution of observed means could theoretically be 
obtained with a variety of hailpad networks of varying qreal 
distribution and given density. 

To illustrate the basics of this problem consider the 


following example (Panofsky and Brier, 1963). When a sample 
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of an apple pie is cut and tasted, it is possible to 
conclude that "This pie is good" or "This pie needs more 
sugar". Of course, neither of these statements is based on 
complete knowledge, since nothing is known about the pie as 
a whole, but only about the sample tasted. What we really 
mean, is, Since the sample chosen tastes good, and since the 
other pieces probably taste the same (i.e. we are making the 
implicit assumption that the pie's properties are 
homogeneously distributed), it is reasonable to presume that 
the rest of the pie will taste good. Thus, from the 
characteristics of the small sample, inferences on the 
Characteristics of the whole are drawn. Actually, one piece 
from the entire pie would be a very poor sample of the pie 
in the statistical sense, for the oven may have been 
Slanted, and all the sugar may have drifted into the sample 
tested. As a resuit drawing conclusions about the whole pie 
would be incorrect. A better method would have been to 
collect a "random" sample by spooning a taste from various 
parts of the pie. In order that all of the sample pieces do 
not come from the same sector of the pie, it is better to 
choose a representative sample from each area. This is 
"stratified sampling"; that is, the pie is divided into M 
uniform sectors and k spoonfuls from each sector are chosen 
at randon. These basic ideas can be applied ina similar 
fashion, using the Monte Carlo method, to impact energy 
density values calculated from hailpad data. Here sampling 


is carried out for impact energy density values on a 
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hailswath basis and not on a hail day basis where more than 
one hailswath may have occurred (i.e. one pie at a time and 
not many pies which may have been cooked on other stoves). 
The method of simulating statistical trials (or 
Sampling) using computational techniques and recording 
numerical characteristics obtained from these experiments is 
called the Monte Carlo method. The term "Monte Carlo" 
originated in 1942, at Los Alamos, in studies conducted by 
von Neumann, and it appeared in publication in the works of 
N. Metropolis and S. Ulam (1949) and Shrieder (1966). For 
impact energy density values the method is to select a 
random sample many times using computational techniques and 
to record the numerical results. The prime requirement of 
the Monte Carlo method is the availability of random 
numbers. These can be obtained on the computer by employing 
one of a number of subroutines which generate random 
numbers, based on arithmetical algorithms from the computer 
library. In this analysis the subroutine used to select the 
samples (i.e. to generate the random numbers) is GGU3 found 
in the *IMSLLIB library. The numbers obtained from this 
subroutine are uniformly distributed on the interval 0,1. 
Figure 3.1 eweee the distribution of 40,000 random numbers 
obtained from subroutine GGU3 for an expanded interval 0,100 
in five unit steps. The results show that the number of 
values in each interval appears to be uniformly distributed 
with a mean of 2000 and a standard deviation of 33.9, with 


the lowest interval having a value of 1918 (4.1 percent 
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below the mean), while the largest interval has a value of 


2066 (3.3 percent above the mean). The above subroutine is 
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Figure 3.1 Distribution of 40,000 random numbers drawn from 
subroutine GGU3 for the expanded interval [0,100] 
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number of ways ! of selecting the random numbers on the 
computer for use in the Monte Carlo method. 

The Monte Carlo experiments for the hailpad data can 
give an estimate of the distribution of the possible average 
impact energy densities which may be observed with a network 
of given size. In theory, if the size of the network (i.e. 
number of stations in a fixed area) becomes infinite, the 
answer to how the Aieeci vat ten of the average impact energy 
densities should behave is given in the central-limit 
theorem of probability theory. This theorem (see Feller, 
1950; pp191ff.) states that the distribution of the sum of k 
mutually independent random variables, with a common 
distribution, approaches a normal distribution as k—+% »« . 
The mean is equal to that of the common distribution and the 
variance equal to k times the variance of the common 
distribution. Hence, for a random sample of impact energy 
density values the distribution of the mean of such a sample 
is expected to approach the normal distribution 2 as the 


sample number increases (ieee. as the number of stations in a 


1 Some of the others ways are: GGBIN-one binomial pseudo 
random deviate; GGBTA-Beta random deviate generator; GGCSS- 
generator of a chi-squared deviate; GGMEXP-generator of 
exponential deviates and GGPOSH-generator of Poisson random 
deviates, among others. 

2 A sample can be tested if it comes from a normal 
distribution if the z score (distribution of the 
standardized variable) of its statistic lies outside the 
range -1.96 to 1.96 for the region of Significance at the 5 
percent level. For the z score distribution the mean is 0 


and the variance is 1. 
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fixed area increases). In practice with the hail energy 
density distributions, it is found that with k greater than 
25, the distribution of the mean impact energy density shows 


some resemblance to the normal distribution. 


3.2 Sampling Techniques 
There are numerous ways of selecting a sample from a 

set of values. In this dissertation five techniques are 

given and described below. Some of these techniques 
represent more the sort of sampling which hailpads do, than 
others. However, all five will be considered at least 
briefly, so that an intercomparison can be made. The five 
technigues are: 

a) The original data are picked at random without 
calculated. Thus, the samples in a single draw are 
Bieeenine a by the elements selected. In this technique 
the number of samples drawn is limited to the number of 
members in the data set or population, since the first 
cannot exceed the latter. 

b) The original data is picked randomly with replacement. 
Here, it is possible to choose the same data value many 
times. This approach has the advantage over the one in 
(a) in that the sample number torr longer restricted 
to the number of members in the data set, but can be 
any size. 


Cc) A distribution is fitted to the original data and 


Teret© Lisd ods Atiw epitoe 1g ae bs 5 3 : 
osdt wersetp A dtiw sate bavo? ai ; $b 4a obs ida: 
 ewode ytkeneb yprsae toeqm asom odd 20 ae haudbase i 
-soltudi ta tb ioateg” s sak ham 


5 iio i 


a 


6 mO1k aliqmse 6 paitosise to oyeH ue 
S18 2oupindoes ‘ovii noisesxeeekb ius Po 2 

_ Beupindses seett to Smoe soled | be iad ixz 
isd3 ,ob ebsytind aoe iw patiquse ro d108 ous “¢ 
tesol t5 borghiremoo od Ebbw pe is . x | 


oviit sd? .obsm od abo noaizsqaooratak te teas o2 atau 
suodsiw mohasi $6 bedoiq S15 arab 
Sis. atstemstsq biieiiaadinds avorzéy | i ag \d sr 
ots wstb sipate 2 nt aolgmse ods aa Be +e 
oupindoes aids aI .bosoeloe eangmoro ed vd bent 
to aedmun oft oF bottmil ai gusab aolqmse Yo. ed 
Ja2xLt edt sonte i ail ht 10 362 BI55 ond ni ae 
F stots l eds ae Wot 7 
-taensvsiqes situ yimobasa Besozg | ai stsb Lsazpiso wae (d 
yas sulpy stsb. ona, edt se0od5. os sidieaoq ak tt — me 
ak 900 out z6v0 opatmrsvhs iy wee dosoxgge aba? ,aomts 
bedoistass zopaoi on ai nome elguse sit tea ai (8) 


aa a aon ee Bieh odgt nr eo jedaon ads oF 


a 


oy ae | ee eee | s9at2 yas 
‘ned ona Hab sentyase eittlos eek SPwobiwdiseth gp 


d) 


€) 


41 


random selection is carried out from the fitted 
distribution. This procedure allows the selection of 
many more random variables than the original data set. 
However, the results are only as reliable as the fit of 
the distribution to the original data. 

impact energy density values from the original data is 
divided into M equal segments of probability values and 
a random selection of k values is made from each 
segment. If M=1, the technigue is equivalent to 
technique (b). If, however, three segments (M=3) are 
used then one-third of the impact energy density values 
will come from the low band of values, another third 
from the intermediate values, while the last third will 
be drawn from the high band of values. This procedure 
guarantees that all the samples do not come from the 
same band of values, as is possible in procedure (b). 
Spatial Saupling: In this techniqueva “true rspatial 
distribution is estimated by an interpolation, 
contouring, or an objective analysis method. This 
technique accounts for the fact that the impact energy 
density values are distributed in a two-dimensional 
space, as is the sampling network. Thus emphasis is 
placed on the spatial distribution of the impact energy 
densities. A random network (but with a certain average 
spacing and hence a fixed density) is imposed and a 


value from each network grid is randomly selected. 
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In Chapter 4 these five techniques for selecting impact 


energy density values are applied to the 29 July 1974 
hailswath as an example of the procedure. The computer 
program for the first four techniques was developed in 
FORTRAN IV on the University of Alberta AMDAHL 470V/6 
computer. It appears in Appendix I. The computer program 
used for the contouring technique in the fifth sampling 
approach was developed by R. Newton (1973) in a modified 


FORTRAN IV language on a DEC PDP-9 computer. 
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CHAPTER 4 
APPLICATION OF THE MONTE CARLO METHOD TO 
HAILPAD DATA 
In this chapter the Monte Carlo method is applied to 
the impact energy density values and the distribution of the 
means of k samples is determined. The method is used with 
the five techniques described in Chapter 3: sampling without 
replacement from the observed distribution of impact energy 
densities; sampling with replacement from the observed 
distribution; sampling from an idealized two-parameter gamma 
distribution; stratified sampling from the observed 
distribution of impact energy densities; and spatial 
sampling from a contoured pattern of the observed energy 
density field. The data from the 29 July 1974 storm are used 
to illustrate the techniques. This was a severe hail day on 
which 34 hailpads were returned in a hailswath area of 3183 
km2. The storm began west of Rocky Mountain House in the 
middle of the afternoon and moved south-east towards Three 
Hills. Figure 4.1 shows the precipitation reports 4 and the 
location of the hailpad observations for this hailswath. 
Golfball-size hail was reported by the volunteer observers 


south of Red Deer. 


Figure 4.2 shows the distribution of the impact energy 


1 In Figure 4.1, the numbers indicate hail of different size 
categories: 1i-shot; 2-pea; 3-grape; 4-walnut; 5S-golfball; 6- 
greater than golfball; 0 is hail of unknown size; R denotes 
rain only and X no rain or hail. 
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density 2 calculated from the 34 hailpads for this 
hailswath. In this figure the impact energy density values 
are grouped in classes of 25 J m. The population mean is 
84.19 a, the standard deviation 120.6 J my the range from 
0.5 to 546.8 J m3 , and the median value 41.2 J m. The 
curve through the triangles in Figure 4.2 is the best 
fitting two-parameter gamma distribution, which is discussed 
in Section 4.3. The cumulative frequency of the energy 
density is displayed in Figure 4.3 where the curve through 
the triangles is again the best fitting two-parameter gamma 
distribution. Other hailswaths were examined and the impact 
energy density distributions for these are presented in 


Appendix III while the results are discussed in Chapter 5. 


4.1 Sampling Without Replacement 

A sample of k impact energy densities is drawn randomly 
N times from the total number Z of impact energy densities 
for the hailswath. The standard deviation of the sample 


mean, denoted by s, is calculated according to: 


2 The impact energy densities for the 29 July 1974 hailswath 


are given in Appendix II. ‘ ; 
3 In three years of data (1974-1976). the hailpad impact 


energy density values ranged from 0.1 to 2332.6 J ms 


2i assem shpat: od? . Tet as. a aseasts ag dequoz « s 
mort opasa edt a &b 8.08? sotisiveb | 7: | ate fe oy if . 
edt Smt Sra oulsy asibea odt heme, Done ilsna iid. 
geod oft 2i $.8 upp? ak aslpasias edt dpuoids ovis 
fsaenoa ib af doidw SEO MAANE BSE sumsy sotomstsq-ows att: 
ypreaes eit to ea aed ov kgahionans edt sae roe : 
dpsords evapo ods ovedw £.8 saupit at boxargut of +iene i 
SMM HP redomsnsq-ows priststit tasd sd? aisps ak asipasizs oat 
gosqmi odt Bas bagimexs ose attswalisd rei20 aseeieaniat:s 
at potaseerq 915 seeds Oz eaoisudiztaLb witenob\reseae 


2 tetqedd al beesvoekb ors ete ad? oLidw : eee xibaoqqa 


| snansontqos tuoditW vatignse T.¥ i 

{imobass dwsrbh, at set tiensb Wasds sougad x zo esiqnse A q 
esitiansh yprsas tosymi to's redeva pee end word eomis u i 
olquse sit Xo apioeeeise busbuste ei? wistswalted sit rok i 


207 pnaibsocos botsltaiso ek ,ea yd fhetoneb 169m 


| 


epee, ay . OTe oes 
(r.b) Cee Ry Tee me 


46 


20 


op] 
ao 


ao 
i=) 


eS 
Oo 


JUL 29 74 


WwW 
Oo 


PERCENTAGE FREQUENCY OF OCCURRENCE 
Nn 
i=) 


_ 
oO 


0 100 200 300 400 500 600 730 
ENERGY DENSITY IN JOGULES/MxM 


Figure 4.2 The percentage frequency of occurrence of impact 

energy density values for 29 July 1974 hailswath (The curve 

through the triangles is the best fitted two-parameter gamma 
distribution) 
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Figure 4.3 The cumulative frequency of impact energy 
density for 29 July 1974 hailswath (The “stars" represent 
the ofiginal data, while the "triangles" the best fitted 

two-paramet@r gamma distribution) 


In Equation 4.1 the sample means are x; (where i=1,Nj) * , 


while the mean of the sample means is given by x. The k 


Samples are drawn at random and the distributior of the means 


She 4 ge oe eee 
* N-1 is used in preference to N in the denominator of 
Equation 4.14.1 in order to obtain an unbiased estimator. 


The difference for the values of N used in this thesis is 
small. 
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of the samples is determined. Using N<<1000 might not 
give as reliable an estimate of the distribution of the 
means while N2771000 would consume too much computer time, 
consequently N=1000 was chosen as a compromise. For the 29 
July 1974 hailswath the total number of impact energy 
densities is Z=34. For sampling without replacement k (the 
Sample size) was chosen as 5, 15 and 25. 

For this hailswath the probability density functions § 
of the sample means are given in Figure 4.4. For k=5, the 
distribution of the sample means varies from 5 to 175 J m% 
with no clearly defined single maximum near the population 
mean of 84.1 J m-. The standard deviation of the sample 
means is 49.3 J m*. For k=25 the distribution becomes 
narrower, with the sample means ranging in value from 45 to 
110 J wm (standard deviation of the means is 12.1 J mw’), and 
it approaches a delta function as k-Z. 

An alternative way of displaying the above information 
is presented in Figure 4.5. Fqui-probability isopileths of 
cumulative frequency of occurrence are drawn as a function 
of the sample size. Figure 4.5 depicts the probability 
distribution function. For k=5, 2.5 percent of the sample 


= -2 ‘ 
means are less than 20 J mor greater than 195 J m. This 


implies that with k=5 there is a 5 percent probability of 


S The impact energy densities are grouped in 5 J m 
intervals (e.g. 0.0-4.9; 5.0-9.9; and so on) in the graphs 

of the probability density function of the sample mean. The 
value plotted is in the middle of this class interval i.e. 

2,5, Iso ands so on. 


2 


foo tipi Oo0r >a pata 
ott Io notsodiwerh oft 20 oad 
,9ait rat equoD) dowm oot sevens 
eS oft 304 se. Lmoxqnoe BRB | | o. 
ypiens sosqat to asditua Lstod eit | 

OAS) 4 snomoontgot tnodtiv pict gape 298 sates ot 4 
tS ae OF ae 28 astiod> 280 fonts 9 

2 entotsoapt ytianeb ysiLidadorg oat lysvelied atdd 
ait, ,€=A TOR oH, i otupid nt wewkp ets anaom sie 
a & Ott oF- = wort 2olisy) atson tenes. ode ao» 
woitsiuvgog S84 1898 abdlixse atpate boalieh ¥ 24 
giqnse ad? to nohI5EVe bosbaste out Pity ren 3 
 Bemigoed xottudt Tekh ang eae 208) Seow e e. 82 2 

o¢ ch mont aphev at ‘palpnss icles Ronet: adie 
his , a & Pash ak io Bias 


foktsnaotal ovods oat on taetaseb 20) yew | 
te edtelgori cae thdwtonysts Be 4 s. Bs 


» 


rotsonvd & ri awsxb § ; th evi 7 oe 

Uithbdddorg als po 395 2. vay sabe Nagle olguze mm: to 
ee 4 
eign se ent to trobaag 2-8 eee, 107 -noxiono? aobsug, ) 


enig Sw eer asa+ “so dseap ‘20 & ps andy ‘pel on 2069 


ae ie 9 wen 
So ‘ean hadozy snsoaeq @ & ei e194s ex date pods eobiqui — 


49 


K 25.0 PGs 
K 216.0 yyy 
70 K 226.0 xx 


n 
io) 


uo 
(=) 


JUL 29 74 


WITHOUT REPLACEMENT 


> 
i=) 


N= 1000 


PERCENTAGE FREQUENCY OF GCCURRENCE 


0 25 50 715 100 125 150 175 
ENERGY IN JGULES/MxM 


Figure 4.4 The probability density function of the sample 
means for 29 July 1974 hailswath (sampling without 
replacement) (Impact energy density values are grouped in 

5 J un classes) 


obtaining a sample mean which is at least 132 percent in 

error from the population mean. AS the sample size increases 
towards the population size of 34 the equi-probability lines 
converge toward the population mean. Figure 4.6 displays the 


standard deviation of the sample means as a function of the 


pees | 
wp ehh 


eT eS JUL ; 
THANIIAIGSA. TUOHT IW , 


= soot =” vi coke hk iin ‘ : ia 


50 


JUL 29 74 


RKITHOUT REPLACEMENT 


pee N= 1000 


26 aaa 
BG) Bac 
40-0 xxx 
50.0 gee 
EES 
TE LOxie x 
1-8 7277 


160 


140 


120 


100 


ENERGY DENSITY 


80 


60 


40 


10° 


Figure 4.5 Isopleths of the probability distribution 
function of the sample means for 29 July 1974 hailswath 
(sampling without replacement) 


sample size. The standard deviation decreases approximately 
as KY’ for increasing values of k from 5 to 25. As the sample 


size increases towards the population size (k=Z=34), the 
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Figure 4.6 The standard deviation of the sample means as a 
function of the sample size (sampling without replacement) 


standard deviation must decrease more rapidly towards zero. 
The dependency of the standard deviation on k can be 
calculated by fitting the data points © with a power-law 
equation of the form 


p_Leseertecpesse, Ff 
6 A least squares fit is applied. 
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standard deviation = (4.2) 
where bis a constant. From Figure 4.6 (for k from 5 to 25) 


this equation gives x=0.87 and b=210. 


4.2 Sampling With Replacement 

Sampling with replacement is similiar to sampling 
without replacement except that once an element has been 
selected it is returned to the population before choosing 
the next element. Hence it is possible to choose the same 
element more than once. The random selection and the 
statistical calculations are carried out in the manner 
described in Section 4.1. The\ sample size is no longer 
limited by the population size (number of hailpads in the 
hailswath) as is the case for sampling without replacement. 
For the storm of 29 July 1974 the graph of the probability 
density function with k equal to 15, 60 and 105 is given in 
Figure 4.7. The distributions in this figure appear to 
approach a Gaussian form as k becomes large. This technique 
is governed by’the central-limit theorem, as compared to 
sampling without replacement (which is not). Figure 4.8 
gives the probability distribution function for the 29 July 
1974 hailswath. In this figure, a sharp decrease is observed 
in the 2.5 and 9755 percent probability isopleths top K=15 
to 60, while the 25 to 75 percent probability isopleths show 
a lesser decrease. The standard deviation of the sample 


means (Figure 4.9) shows a steady decrease with increasing 
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sample size (for k=15 to 105, x=e51 and b=125 ineequation 
4.2). This is in good agreement with the central-limit 
theorem, which gives x=.50 and b=120.6. Im Figure 4.9 is 
also shown the variation of the standard deviation of the 


means as a function of the sample size for k=1 to 15. For 


athed 


Av eS JUL : an 
THIWASRISSA HETK 
2001 


=u 


ro 


y. 
a. 


200 ; : 
JUL 29 74 2.3. daa 
26-6 444 
40.0 xxx 
WITH REPLACEMENT 60-0 »66 
80.9 a46@ 
76.0 ERX 
180 N= 1006 O-6r2Z 
160 
140 
120 


— 

— 

tA 

w 

ea 

iJ x 

aa é 

cD 

a ———_——_—__ 

ud aaa 

az ee 

. : ee 
ie Bar # 
40 


20 


NUMBER 


Figure 4.8 Isopleths of the probability distribution 
function of the sample means for 29 July 1974 hailswath 
(sampling with replacement) 
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Figure 4.9 The standard deviation of the means as a function 
of the sample number (sampling with replacement) 


k=1, the standard deviation equals 126.5 J n’, * (use 
Slightly higher than 125 J m-~ the value computed for k=15 to 


105 using Equation 4.2, and higher than 120.6 J n~* (the 


7 The standard deviation equals 126.5 J m~* because N=1000, 
as N-»yooit is expected that this standard deviation would be 
equal to _the standard deviation of the population mean. 
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standard deviation of the population). 


4.3 Sampling From a Fitted Two-parameter Gamma Distribution 

The third procedure is to fit a two-parameter gamma 
distribution to the original data and apply the Monte Carlo 
method by sampling from the fitted distribution. Other 
distributions can be used, if a good data fit can be 
obtained (e.g. Kappa 2 distribution). Before the results 
from the data of the 29 July 1974 are presented, the gamma 
distribution is discussed. The impact energy density values 
derived from three seasons of hatlpad data are positive and 
range in value from Obi x 2332.6 J wm. Low impact energy 
densities occur most frequently with a rapid decrease to 
high impact energy density where the relative frequency of 
occurrence is small. A distribution of this nature can be 
fitted quite well using the two-parameter gamma distribution 
(see Figure 4.10). This distribution has been found useful 
in fitting and transforming distributions of various 
meteorological variates which are restricted to positive 
values such as precipitation, vapor pressure and 
precipitable water (Panofsky and Brier, 1963). 


The two-parameter gamma distribution (Keeping,1962) is 


defined by the frequency function: 


(4.3) 
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Figure 4.10 Density function for the gamma distribution 


where /7( 2) is the gamma function defined as: 


P(x) = Gs (4-4) 
° (x 70) 


and 


Ma) = (c -1) 1 («-1) = j-)4- 5) 


if «< is an integer. 
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The mean of the gamma distribution is XB and its standard 
deviation is pie ° 
The optimal estimates of x and pe (Panofsky and 


Brier, 1963) using: 


and 


epi @ 
p= = (4.8) 


Using Equations 4.7 and 4.8, the two parameters ({ ~ and 
A) for the population data 8 were calculated and a 
frequency distribution obtained. In this technique it is 
possible to select from the distribution an impact energy 
density value which does not exist in the original data 
base. In Figure 4.2, the fitted two-parameter gamma 
distribution is given by the curved line, with triangles in 
the middle of the class interval. The fit of the two-- 
parameter gamma distribution as a cumulative frequency is 
demonstrated in Figure 4.3. Figure 4.2 shows that for impact 
energy values less than about 100 J nthe gamma 


distribution slightly underestimates the frequency of 


. 7 " 
® For the 29 July 1974 data: % =.63, f =133.5, af (mean) =84.1 
and Al (standard deviation) =106.0. 
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occurrence of the original data and why the standard 
deviation obtained from the gamma distribution is less than 
that of the population. Grouping the impact energy densities 
FRM cs n-* classes for the 29 July 1974 hailswath the value 
of X was 6.17. With V5¢=21-0 (with 12 degrees of freedom) 
we conclude that the gamma distribution gives a good fit to 
the original data. The mean of the two-parameter gamma 
distribution is 84.13 J m” * (same as the population mean) and 
the standard deviation is 106.0 J m * (compared with the 
population mean of 120.6). 

The Monte Carlo method is applied using k equal to 15, 
60 and 105. Figure 4.11 tose the distribution of the sample 
means for the various sample sizes. For k greater than 15 
the distribution begins to exhibit a Gaussian form as 
expected according to the central limit theorem. Unlike the 
technique of sampling without replacement, in this procedure 
k is no longer restricted by the size of the population. The 
disadvantage of this procedure is that it requires a good 
fit to the original data, while for the other two procedures 
no such requirement is necessary. Comparing the probability 
of frequency of occurrence of impact energy values for given 
sample numbers (Figures 4.5, 4.8 and 4.12), it is observed 
that for the fitted two-parameter gamma distribution the 
sample means will be within 25 percent of the population 
mean more frequently than in either of the previous two 
procedures. The 25 and 75 percent probability isopieths are 
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Figure 4.11 The probability density function of the sample 
means for the 29 July 1974 hailswath (sampling from a two- 
parameter gamma distribution) 


similarly observed for sampling with replacement. Figure 
4.13 ssh the standard deviation of the sample meass as a 
function of sample size. The standard deviation of the 
sample means shows a steady decrease with increasing sample 


size (for k=15 to 105, x=.54 and b=120 in Equation 4,2) and 
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Figure 4.12 Isopleths of the probability distribution 
function of the sample means for 29 July 1974 hailswath 
(sampling from a two-parameter gamma distribution) 


shows agreement with the central-limit theorem and the work 


of Strong and Lozowski (1977). 


4.4 Stratified Sampling 


In stratified sampling, the cumulative distribution of 
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Figure 4.13 The standard deviation of sample means as a 
function of the sample size (sampling from a two-parameter 
gamma distribution) 


the impact energy density is divided into M equal 
probability segments. From each segment k/M 9 samples are 
randomly selected with replacement using subroutine GGU3. 


9 k/M should be an integer, to assure that the same number 
of samples are obtained from each of the M segments. 
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“the total of k samples are used to give the sample mean. The 
process was repeated 1000 times for the 29 July 1974 


hailswath. This technique has the advantage over the 
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Figure 4.14 The probability density function of the sample 
means for 29 July 1974 hailswath (stratified sampling) 


previous three techniques in that the sampling is made from 


various impact energy density bands. In the other approaches 
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(it is possible that all k values could come from either a 
low or high energy region. With stratified sampling this is 
not possible. It is therefore believed that this type of 
Sampling may better simulate the kind of sampling performed 
by the hailpad network for those hailswaths where the impact 
energy distribution is not patchy. 

For the example of the 29 July 1974 storm, two cases 
are presented: (1) M=5 (five segments) with k (the total 
sample size) being 15, 60 and 105 and (2) M=15 (fifteen 
segments) with k having the same values as in the first 
case. The graph of the probability density function for M=5 
is Figure 4.14 while for 4=15 it is Figure 4.15, The 
distributions are narrower than for non-stratified sampling, 
and show a resemblance to the Gaussian form. The 
distributions for the same k become narrower for M=15 
(Figure 4.15) than is observed for M=5 (Figure 4.14). Figure 
4.16 shows that 95 percent of the mean impact energy 
densities lie within a range of 42 J Ww (for M=5) as compared 
to within 13 J m* for M=15 (Figure 4.17). The equi- 
probability lines are clustered in a very narrow impact 
energy band in Figure 4.17 as compared to that observed in 
Figures 4.8, 4.12 or 4.16. The standard deviation =f the 
sample means for M=3, 5 and 15 are shown in Gio ire GeMS In 
all tires cases, the standard deviation of the saiuple means 
decreases as k increase from 15 to 150. The values 9f x for 
M=3, 5 and 15 are .50, .49 and .49 (in Equation 4.2), while 


the values of b are 62.0, 42.0 and 14.7 respectively, 
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Figure 4.15 The probability density function of the sample 
means for the 29 July 1974 hailswath (stratified sampling) 


4.5 Spatial Sampling 

In this section, the attempt is to more closely 
approach the network sampling by trying to simulate two- 
dimensional spatial sampling computationally. In this 
technique a "true" spatial distribution is estimated by a 
contoured interpolation method. This technique takes into 
account the fact that the impact energy density values are 


distributed in two-dimensional space, as is the sampling 
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Figure 4.16 The probability distribution function of the 
sample means for 29 July 1974 hailswath (stratified 
sampling) 
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Figure 4.17 The probability distribution function of the 
sample means for 29 July 1974 hailswath (stratified 
sampling) 
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Figure 4.18 The standard deviation of the sample means as a 
function of the sample size (stratified sampling) 


network 19 , The contouring programs used in this 


dissertation are those developed for a DEC PDP-9 computer by 


10 In reality, the hailpad sampling network is distributed 
in three dimensions, because of a variation in height above 
sea level from station to station. Here it is assumed that 
this variation is minimal. This is a valid assumption for 
central Alberta. 
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R. Newton(1973) 12 , 

The first approach is to estimate the "true" 
distribution by using a method of statistical prediction 12 
Res kaka by Newton(1973) for the 34 impact energy densities 
of the 29 July 1974 hailswath. The results using the 
contouring program are shown in Figure 4.19 13 , 

Next in order to simulate two-dimensional spatial 
sampling computationally a square grid is imposed on the 
two-dimensional contoured impact energy density analysis and 
points are chosen at random within each grid square. The 
first case considered was a square grid of length 11.5 km 
(132 km2) !* and from each grid square {k=77) through which 
the hailswath crossed, an impact energy density was chosen 
at random !5 . The impact energy densities were then added 
ae the grids, and the simulation repeated one thousand 
times. Next, a larger area was selected such that the length 
of the square grid was 22.9 km (526 km@). Again from each 
grid through which the hailswath crossed, an impact energy 


aoe 


11 Unfortunately its conversion to the University of Alberta 
AMDAHL 470V/6 is still pending. Hence, the examples given in 
this section are the result of extensive manipulation and 
data transfers between the two computers. 

12 A technique in which a "probability function" is used to 
predict the "most probable" value at any point on a map area 
which:has been defined by a number of sampled data values. 
From a practical point of view this might be considered a 
system of data averaging. 

13 In Section 4.5, the north direction on the graphs is up, 
while west is to the left-hand side. 

14 Figure 4.20 shows a square grid of length 11.5 km 
superimposed on this figure. 

1S An x-value was first chosen at random using GGU3, then a 
y-value was similarly chosen. 
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Figure 4.19 The contoured values for the 34 impact energy 
densities from the 29 July 1974 hailswath. A square grid of 
length 11.5 km is superimposed on the contoured swath. The 
contouring is based on actual values at the hailpad sites, 

however, for the purpose of displaying the impact energy 

densities in graph form the contoured values for this 
hailswath were assigned unit values from 1 to 9 with 1 
representing the 50.0 J m“*contour ; 2 the 100.0 J 
m-* contour; and so on. 


density was chosen at random. For the 526 km? grids, k was 
14, choosing one impact energy density in the grids 
influenced by the hailswath. Figure 4.20 shows the 


probability density function obtained using the 132 km? and 


ee 
of i 2 ies 
\@ et ae a oe 


ame aes a i ‘oe : BP ie Maas ine) ne f 
; ye, yee ate: a> i se 


f 


ut eae, fy f 

RE let 
+e Hy e's ele ; Midi pe «4 ; a My 
: be i 


i 
+ 
reg 


ee 
A. hd ee 


a ci a) 
Sse AA 
asw a) 

a r = ® be 
7 q 


71 


70 


on 
i=) 


2a] 
i=) 


JUL 29 74 


2-0 CONTOUR 


pS 
oO 


N= 1600 


PERCENTAGE FREQUENCY OF OCCURRENCE 


0 25 50 75 100 125 159 175 
ENERGY IN JOULES/MxM 


Figure 4.20 The probability density function (spatial 
sampling) 


526 km2 square grids. The distribution shows a resemblance 
to the Gaussian form with an increase in sample size. The 
standard deviation of the sample means as a function of the 
sample size (given in Figure 4.21) decreases with increasing 
sample size (x=.36 and b=32 in Equation 4.2). The slope of 


the standard deviatiam is lower than that obtained for 
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stratified sampling (for M=3,x=.50) as is also b (the 


standard deviation intercept). 
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Figure 4.21 The standard deviation of the sample means as 
function of the sample size (spatial sampling) 


In summary, of the five procedures, sampling without 
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replacement is the only one in which k «the sample size) is 
restricted to the population size. It is questionable 
whether the procedure itself represents the sort of sampling 
that can be expected from a hailpad network due to the 
restricted sample size. The second technique, sampling with 
replacement, uses the original data which seems to be 
satisfactory; however, this method suffers ae the 
disadvantage that the k values may come from either a low or 
a high band of impact energy densities. If such is the case, 
a bias can occur in the extremes. In order to overcome this 
disadvantage, stratified sampling is tested as the fourth 
procedure. Both cecudaltas show the same decrease in slope 
in the standard deviation of the sample means, however, for 
stratified sampling the standard deviation intercept (b in 
Equation 4.2) is lower than that for non-stratified 
sampling. Also for stratified sampling, the sample means 
show less variation in impact energy density (in the 
probability distribution function) than that observed for 
non-stratified sampling. The third technique introduced is 
that of sampling from a fitted distribution. This approach 
is satisfactory if a good fit to the original data is 
obtained. The results 16 from this technique are comparable 
with those obtained for sampling with replacement (non- 


stratified). The last approach is spatial sampling, and 


16 The probability density function, probability 
distribution function, and standard deviation for a given 
sample size. 
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since the areal distribution of impact energy densities is 
considered in this technique, it is believed that this 
approach best simulates the type of sampling performed by a 
hailpad network. In all the techniques examined (except 
sampling without replacement) the standard deviation of the 


Yo = 2 
« The relationship 


sample means behaved approximately as k 
for sampling without replacement was kK "’ for k from 5 to 25 
with the exponent of k increasing rapidly to a liarge value 
17 as k approaches 34 (the population mean). 

In this chapter five procedures were discussed using 
the data from the 29 July 1974 hailswath. Other hailswaths 


were analyzed in a similiar fashion and the resuits for 


these are given in Chapter 5. 


17 The slope approaches ~~ as k approaches 34. 
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CHAPTER 5 
RESULTS FOR OTHER STORMS 


5.1 Introduction 

In this chapter, the results from five hailswaths 4 in 
the three years 1974(3) 2 and 1975(2) are examined using the 
first four techniques described in Chapters 3 and 4. At 
present the fifth technique requires that the DEC PDP-9 
computer programs (developed by Newton) be transferred and 
nodified to the University of Alberta AMDAHL 470V/6 
computer. Further extensive computer program development is 
required as support for the modified computer programs 
before the technique could be used for analysis of hailpad 
data. This technique, it is believed is the sort of sampling 
that hailpads do in the network and therefore the technique 
is suggested as an area of future research. Only impact 
energy densities within the hailswath are used in the 
analysis. In most cases other impact energy densities are 
available for the hail day, but some of these are not 


associated with the hailswath 3 since they were produced by 


1 The hailswath of the 29 July 1974 is included as the fifth 
hailswath in this dissertation. 

2 The number in parenthesis is the number of hailswaths used 
for each year. 

3 Precipitation and hailfall reports and radar echo 
information for the hail day were used to determine which 
impact energy density values were associated with the 
hailswath of interest. 
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other storms occurring on that day. Since hailswaths have 
different characteristics (for example, motion and direction 
of the hailswath, maximum hail size they produce, among 
others) it is not unreasonable that impact energy density 
distributions also vary among hailswaths. Hence impact 
energy densities from the other storms which cccurred on the 
same day are not used in the analysis. 

The maps of hail and precipitation observations + 
reported by the volunteer observers are found in Appendix 
II. Six graphs resulting from the application of the various 
techniques to the analyzed hailswaths are shown in the 
appendices 5% . The graphs of percentage frequency of 
occurrence of impact energy densities for each hailswath 
analyzed are found in Appendix III, while those for the 
cumulative frequency are in Appendix IV. The other four 
graphs (probability density function; probability 
distribution function; standard deviation of the sample 
means as a function of the sample size and standard 
deviation of the sample means versus the number of hailpads 
per 1000 km2) are found in each of the following four 
appendices: 

a) Appendix V: Graphs for Sampling Without Replacement 


* These maps were plotted on the Calcomp Plotter connected 
to the University of Alberta computer by a computer progran 
(unpublished) developed by the author in 1968 and modified 
in 1974 and 1975. 

S Results for the 29 July 1974 hailswath are not given in 
the appendices since they appear in Chapter 4. 
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b) Appendix VI: Graphs for Sampling With Replacement 
Cc) Appendix VII: Graphs for Sampling From a Fitted Two- 
parameter Gamma Distribution 
d) Appendix VIII: Graphs for Stratified Sampling 
Before we examine the four procedures, a summary of the 


hailswaths © analyzed is given in Table 5.1. The first 


Date No. of Area Max.Size BV (Sci s OL ak 
Pads ({km2@) cm Density Density 
24 June 1974 20 3820 SU 107.8 1412.0 
5 July 1974 21 1303 1.4 81.6 70.6 
29 July 1974 34 3183 220 84.1 120.6 
Selubypet?75 18 886 Bs 131.6 124.0 
20 July A975 32 883 Zo TOS 101.0 


column gives the date of the hailswath while the second 
column indicates the number of dented hailpads used in this 
study. The third column gives the area (km?) of the 
hailswath in which the hailpads were collected; this area 
was determined from hailpad and hailcard observations. The 
boundary of the hailswath was sketched by hand and the area 
calculated using a planimeter. The fourth cclumn gives the 
maximum hail size (in cm) observed by the hailpads, the 


fifth column the average impact energy density, as 


6 The hailswaths considered in this dissertation are 
moderate to large (in terms of area and impact energy 
densities) and hence the generality of the results may be 
limited. 
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calculated from the hailpads while the last column gives the 
standard deviation of the calculated impact energy 
densities. The average impact energy densities given here 
show the typical values that have been found in the hailpad 
data of 1974 to 1976 when there were more than 15 hailpads 
in a hailswath. These results then are representative for 


the larger, more damaging hailstorms. 


5.2 Comparison of the Sampling Techniques 

In this section three techniques ? are examined using 
the five hailswaths. The technique of sampling without 
replacement it is believed does not represent the sort of 
sampling that can be expected from a hailpad network, and 
hence this technique is discussed briefly. In sampling from 
a fitted distribution it is desirable that a good fit be 
obtained to the original data. Consequently we now examine 
the results of fitting the two-parameter gamma distribution 
to data from the five hailswaths (see Appendix III for the 
cumulative frequency graphs). A measure of the discrepancy 
existing between the observed and expected frequencies is 


supplied by the statistic Se (Spiegel, 1961) given by 


‘ gle = $i) (5.1) 


7 The three techniques are: sampling with replacement (non- 
stratified, saupling. trom a fitted distribution and 
stratified sanpling. 
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where Oo; and @; represent respectively the observed and 
expected frequencies in the jth class interval. If x? =0, 
the observed and the fitted frequencies agree exactly; while 
the larger the value of XSF « the greater is the discrepancy 
between the observed and expected frequencies. Here the 
expected frequencies are computed on the basis of the 
bypothéesisthH,*7 that.is if por given by Equation 5.1 is 
greater than XA<(eritical value atthe 05 significance 
level), then we would conclude that the observed frequencies 
differ significantly from the expected frequencies and would 
reject H, at the 5 percent level of significance. Otherwise 
the hypothesis would not be rejected. The results for the 


five hailswaths are given in Table 5.2. The first column 


Table 5.2 Summary of the Chi-Square Test for the Hailswaths 
Analyzed 
= 2 2 ; : : 
Date x asks Kes X99 Description of Fit 


24 June 1974 6.9 "12 721.0 42652 “Good 

Sr uly 1974 18.8 8 15.5 20.1 Reasonable 
29 July 19-74-—662 12 21.0 °26.2 Good 

3 July 197571323 212 G2490e42642 Good 
20 July ADFSCV4EST GI240727.0ha26024 Goad 


gives the date of the hailswath; the second, the statistic 


Xe using Equation 5.1 with the impact energy densities 
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grouped in 25 J m “classes 8 ; the third, the number of 
degrees of freedom; the fourth, the .95 level of 
significance; the fifth, the .99 level of significance; and 
the last column, a description of the fit 9 . Overall, for 
the five hailswaths examined in this dissertation, the two- 
parameter gamma distribution gives a reasonable fit to the 


hailpad data 1° . 


5.2.1 Probability Density Function 

The probability density functions show resemblance to 
the Gaussian form with increase in the sample size (i.e. for 
values greater than 15) for the techniques examined. The 
resemblance to this form becomes more rapid (and hence more 
apparent) for stratified sampling than for either non- 
stratified sampling or sampling from the fitted 
distribution. It is also for stratified sampling that the 


probability density function shows the largest percentage of 


6 The class intervals are grouped such that the lowest 
observed frequency is greater than one unit. "In practice, 
the classes with low expected frequency usuaily come near 
the ends of the distribution, and it is common practice to 
combine or "pool" the end-classes until the ¢; (expected 
frequency in the ith class) reach a satisfactory size. The 
objection to pooling is that some important differences 
between ee and e- in the end-classes may be hidden by this 
treatment." (Keeping, 1962, p.252) zs 

° A fit is described as good if AN< Aes, 3 ceasonable if 
Nie gk OX; Acandipoor tkfi Axi 2905 ° 

10 The two-parameter gamma aistribution implies that there 
is a very small but finite (non-zero) probability of 
obtaining a very large impact energy density (i.e. there is 
no maximum impact energy density cut-off). Thus, in a sense, 
this distribution can be considered as a poor fit to any 
hailpad impact energy density distribution. 
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occurrence for a given sample size. 


5.2.2 Probability Distribution Function 

The graphs of the probability distribution function are 
another way of displaying the results found in the graphs of 
the probability density function. For a given hailswath, the 
spread in range of the 95 percent confidence limits varies 
from technique to technique. The equi-probability isopleth 
lines converge towards the population mean with increasing 
sample size for all the hailswaths examined. The greater the 
spread between these lines the greater the variability in 
the impact energy densities. Sampling with replacement and 
sampling from the fitted distribution show comparable 
variability in the 95 percent confidence limits of the 
impact energy density. The improvement in the variability of 
the 95 percent confidence limits (compared to the above two 
technigues) shows up in stratified sampling. Here (for k 
from 15 to 105), the 95 percent confidence limits vary from 
20 to 50 percent of the impact energy density obtained for 


the other two techniques. 


5.2.3 Standard Deviation of the Sample Means as a Function 
of the Sample Size 

For the hailswaths analyzed the standard deviation 
decreases with increasing sample size. It is this rate of 
decrease which is of interest. Using a least squares fit in 


Equation 4.2 it is possible to obtain an estimate of how the 
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Standard deviation of the sample means varies with the 
sample size. The results are shown in Table 5.3. For the 
hailswaths analyzed, the standard deviation of sample means 
varies from x "to Seen for sampling without replacement. 
For the other three techniques, little variation is found in 
the exponent of the power law approximation of the standard 
deviation of the sample means. The variation in the x (in 
oy for 7 in ithe range 15:to 105 is from 2.48 to .54 for the 
three techniques. This shows good agreement with x=.50 given 
by the central-limit theorem. Hence for the hailswaths 


/ 
examined the standard deviation varies essentially as Kk, 


fable 5.3 Summary of exponent in a power law approximation 
to the variation of the standard deviation of the sample 
means as a function of the sample size for the four 
techniques 
Date Without With Fitted Stratified 
24 June 1974 1.06 oot 249 » 48 
Saduiy (419789 102 re, 051 .48 
So July, 1976. 5.87 me 254 249 
aiguLy | 4975 | 1.134 249 249 «50 
20 guiy 1975.  .93 49 -48 -50 


In the above discussion, the sample number was not 
normalized for the area of the hailswaths from which the 
data was collected. Of the hailswaths considered, some 
hailswaths are at least three times as big as other 
hailswaths in area, while the number of hailpad observations 
is only one and half times as great. To compensate for the 


area difference, the sample size is expressed in terms of 
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the number of hailpad observations collected per 1000 km2. 
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Figure 5.1 The mean difference between estimates (in 
percent) for sample density in a 1000 km? area. The solid 


Tine as the best fit, line. 


(sampling with replacement) 
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The graphs of the variation of the standard deviation with 


the normalized sample density are given in the appendices. 


The slope of the standard deviation for the normalized 
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sample density is the same as for the non-normalized sample 


density (values are the same as that given in Table 5.3). 
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Figure 5.2 The mean difference between estimates (in percent) 
for a sample density of 1000 km? (sampling from the best 
fitted two-parameter gamma distribution) 


5.2.4 Mean Difference Between Estimates for a Sampling 
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STRATIFIED SAMPLING 
N= 1000 Bs 


t= 5 


MEAN DIFFERENCE BETWEEN ESTIMATES 


| 2 3 4 © 67885 2 3 4 6 @€@788 
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Figure 5.3 The mean difference between estimates {in 
percent) for a sample density of 1000 km2 (stratified 
Sampling) (M=5) 


Density per 1000 km2 


Figure 5.1 shows the variation of the mean absolute 
difference between the sample means ace (ie) and the mean 
of the sample means x, (in percent) as a function of a sample 
size per 1000 km2 (for sampling with replacement) for the 


five hailswaths examined. This graph is the counterpart for 
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hailfall (for sampling with replacement) to the work of 
Collier et.al.(1975) on rainfall (see Chapter 1, Figure 
1.7). In Figure 5.1 the mean difference decreases with an 
increasing density of hailpads in the 1000 km? area. The 
least squares line (solid line in Figure 5.1) gives x equal 
to .50 (where the mean difference is proportional to kK”) 
for the five hailswaths. Thus, on the average for the 
hailswaths analyzed, a factor of 2 decrease in the mean 
difference requires an increase of 4 fold in the sampling 
density. Similarly, the graph for sampling froma fitted 
two-parameter distribution is shown in Figure 5.2. The value 
O©fn sr ase<s50r £orfthis eile disc also. For stratified sampling 
(with M=5) the least squares curve is given in Figure 5.3 
for the five hailswaths (with x=.50). The slopes of the 
least square lines are about the same in all three 
techniques although the intercept on the vertical axis (mean 
difference between estimates or value of b in Equation 4.2) 
is different. Of the three techniques, stratified sampling 
has the lowest mean difference between estimates, while the 
other two techniques are about equal. 

In summary, to observe the same decrease in the mean 
difference, all three techniques require the same increase 
in sampling density. However, for a specific mean difference 
between estimates, stratified sampling requires a lower 


sample density than the other two techniques. 
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CHAPTER 6 


CONCLUSIONS 

In this dissertation five simulation techniques have 
been presented and applied using the Monte Carlo method in 
an attempt to simulate sampling by a hailpad network. The 
first four techniques have been examined for five hailswaths 
that occurred in the last three years (1974 to 1976) in 
central Alberta. Which technique best simulates the sort of 
Sampling that hailpads do in a hailswath, depends on the 
nature of the hailswath itself. If the impact energy 
densities are of a continuous nature ! within the hailswath 
then the technique of stratified sampling (with replacement) 
might best simulate the sampling that hailpads would do in 
that hailswath. If on the other hand, the number of these 
"pockets" is large compared with the number of hailpad 
stations then a technique like non-stratified sampling (with 
replacement) might best describe the procedure. This then 
brings up the question of when stratified sampling may give 
a better representation of network sampling and when it may 
not. The answer to this question may lie in the auto- 
correlation function for the impact energy densities of the 
hailswath. An integral scale A can be defined by 
1 By continous nature we mean pockets of impact energy 
densities decrease continuously from areas of maximum impact 


energy density to areas of minimum impact energy density and 
then increase to areas of maximum impact energy density. 
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Ls dof obtafiel, Naiety (6.1) 


A~- yr-30 


where cC(x,y) is the normalized correlation coefficient. This 


coefficient is defined as 
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where Eluy) is the energy density field. 

If the mean hailpad separation is much smaller than & , 
stratified sampling sweuuie be a close approximation to the 
network sampling. If on the other hand, the mean hailpad 
separation is much greater than & , then random sampling 
(with replacement) may be as good a representation of 
network sampling. When the mean hailpad separation is about 
equal to & , then it is questionable as to which approach 
(if any) is appropriate for simulating the sampling of the 
hailpad network. 

Sampling with replacement (non-stratified), is probably 
best applicable to hailswaths where the impact energy 
density distribution in any small sub-section of the 
hailswath is essentially the same as that in the entire 
hailswath. Sampling from a fitted two-parameter gamma 
distribution is satisfactory (compared to random sampling 


from the original distribution) if a good fit to the 
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Original data can be obtained 2. The results from the five 
hailswaths show that this distribution gives a reasonable 
fit to the impact energy densities. This technique is 
probably best applicable to severe hailswaths where a large 
range 3 of impact energy densities is observed. 

The results of this study show that the probability 
density function of the sample means approximates the 
Gaussian form * in the four techniques (in most cases) for 
sample sizes greater than 15. The variation of the standard 
deviation as a function of the sample size is on the average 
proportional to k ~ (k in the range 15 to 105) for the 
three techniques S$ . In all three of these techniques the 
standard deviation showed agreement with the central-limit 
theorem, that the sample size is proportional to ae 

The graphs of the variation of the mean difference 
between sample means as a function of the number of hailpads 
per 1000km2 for the five hailswaths are the counter-part to 
those given by Collier et.al.(1975) for rainfall. This is 
the first time that results for hailswaths have been given 


in this type of display. The results show that the mean 


2 Stratified sampling can also be performed for the fitted 
distribution. 

3 Impact energy densities have been observed to vary from .1 
to over 2000 J m*. 

* A quantitative estimate of how well the probability 
density function of the sample means approximates the 
Gaussian form can be obtained by using X . 

5S The three techniques are: sampling with replacement (non- 
stratified), sampling from a fitted distribution and 
stratified sampling. 
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difference between the sample mean and the population mean 
is on the average proportional to k~ “ for the three 
techniques. Hence for the hailswaths analyzed, a factor of 
two decrease in the mean difference requires an increase of 
four-fold in the sampling density. 

The results further show that in order to commit a 10 
percent error © between the sample means and the mean of the 
sample means, 90 hailpads are required in a 1000 km? area 
(or 1 hailpad per 11.1 km2) if non-stratified sampling is 
representative of what the network does. The number of 
hailpads (again for a 10 percent error in the impact energy 
density) increases to 200 per 1000 km2@ area (1 hailpad per 
5.0 km?) using the fitted two-parameter gamma distribution. 
An improvement in the sampling density is observed for 
stratified sampling (M=5), where for a 10 percent error 
between the estimates, only 10 hailpads are needed in a 1000 
km2 area (1 hailipad per 100 km?). 

These results show that the dense hailpad network (see 
Chapter 2) has a sufficient number of hailpads in the 
network to determine the mean energy density of the hailfall 
within a 10 percent error ” . The precipitation network, on 
the other hand ® , has an insufficient number of stations in 


the network to obtain in sampling a 10 percent error in the 


6 The difference between seeded and unseeded storms, which 
one is trying to discover may be as small as 10 percent. 

7 Reguiring a 100 percent response from the hailpad network 
operators in the area. 

8 With 500 stations over an area of 34,467 km2 
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impact energy density % . The minimum error in determining 
the total impact energy density for hailswaths of the kind 
considered here with the present precipitation network would 
be about 18 percent, !° assuming a 100 percent response 11! 


in the hailpad observation network. 


9 For sampling with replacement and sampling from a fitted 
two-parameter distribution. 

10 Using results from the fitted two-parameter gamma 
distribution. 

11 A conservative estimate of the response on a given 
hailstorm would be 20 percent. The non-response is due 
hailpads not having been set out as a result of holidays, 
lack of time and so on. This is not expected to lead toa 
bias in the returns. 
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TABLE OF SYMBOLS 


-alpha (parameter used in the gamma distribution) 


~7number of hailstones in the ith size category 


-beta(scale parameter used in the gamma distribution) 
-drag coefficient of a hailstone in free-fall 
-diameter of a hailstone 

-dent diameter 

-hail diameter 

-diameter of hailstones in the ith size category 


2 at 53° W) 


-acceleration due to gravity (981.3 cm s 
-impact energy 


-impact energy density (without a wind) 


yimpact energy density (with a horizontal wind V,) 


-sample size 

-mass of a hailstone 

-number of segments 

-number of times a sampling technique is repeated 
-density of a haiistone 

-density of air 

-standard deviation of the sample means 
-horizontal wind speed 

-terminal velocity of a hailstone 

-angle to the vertical at which the hailstone strikes 
-sample means 


-mean of the sample means 


Z -total number of hailpads in a hailswath 
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APPENDIX I 


FORTRAN IV COMPUTER PROGRAMS 

The following pages contain a listing of the FORTRAN IV 
computer programs used on the University of Alberta AMDAHL 
470V/6 computer, in this dissertation. Program PAD1 was used 
for the first four techniques in the Monte Carlo method. The 
input values required for initiating the computer program 
are: 

(a) N-number of draws(integer value from 1 to 99,999) 

(b) kK, >maximum sample size(integer value from 1 to 999) 

(c) k-minimum sample size(integer value from 1 to 999) 

(dj) 4S k-increments of k(integer value from 1 to 998) 

(e) A-area of hailswath in km(integer value from 1 to 
99;999) 

(£) Type of technique:01-Sampling from a fitted two- 
parameter gamma distribution; 02-Sampling without 
replacement;03-Sampling with replacement; 04-Stratified 
sampling. 

Program PAD2 was modified from PAD1 to accomplish 
spatial sampling. Input parameters (a) and (e) are required 
for this program. This computer program requires; data which 
has been contoured in a 2-D spaceqhere computer programs 
developed by R.Newton were used and applied on ihe DEC PDP-9 
computer). 

Program PLOT was used to produce maps of tue hailfall 


and the hailpad locations. Here hailpad and hai’fall data 
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are required for initiation of the computer progran. 
All three computer programs given here require support 
of subroutines from various libraries (i.e. *PLOTLIB, 


*IMSLLIB). 
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APPLICATION OF HONTE CARLO METHOD TO 
HAILPAD TMNPACT ENERGY VALUES 

DIMENSION Bea ee eee a) 
DINENSTOWIWAG20)2,531 (120,20), V8 (29) ,CE (9, 20) 

DIMENSION VY¥ (500) ,VX (60) ,V¥ (00) ,CE (500 0) 

DIMENSION. AVT(2C) ,AVB(2C) ,-VRA(60) ,V2B (50) 

DIMENSION PX (60) ,P9(00),PP{(60) , O00 VW) 
DIMENSION PPP(12C),PQP (120) ,IPMM (2) 

DIMENSION VIi (120) ,VJJ (120), VKK (120) 

DINENSION VXXXX(50) ,W3(2000) ,VYYYY (50) 
DIMENSION WwWW(5000) 

DIMENSION XZ2Z(20C0) 

COASNON BUF (2048) 

CALL PLOTS (BUF,8192) 

DATA QQ/**!/ 


TRIS PFOGRAM APPLIES THE MONTE CARLO METHOD USING 
THRE EADIFFERENT SANPLING! TECHNEQUES?S ¥ IE. 1) CODE 001 
SANPELNGEEROMS AS FITTED RHOPRARANETER GAHSA DISTRIB *N 
2) CODE, OO2"SAHPLING WITHOUT REPLACENENT. 


3) CODE 003 SAMPLING WITH REPLACEMENT. 
THE FOLLOWING INPUT PARAMETERS ABE REQUIRED: 
NTT(I5) FORMAT) -IS THE NUGBER OF DRAWS (TYPICAL VALUE 


IS 1000) . 

IWP(I3 FORMAT) -IS THE MAXIUM SAMPLE NUMBER STARTING 
AT 5 AND INCREASING IN INCREHENTS OF 25 

KKK (I3 FORMAT) IS ONE OF THREE SAMPLING TECHNIQUES 
DATA INPUT IS UF THE FORM:A(I, De coe 
3X,¥"8.2,A3,A2,43}) GN UNIT 8 WHERE A(T,3} ARE THE 
IMPACT ENERGY VALUES. AS1,AS2,AS3 IS THE DATE OF THE 
HAILPAD VALUE. OUTPUT REQUIRES THREE UNITS: 

9 UNITS FOR THE PLUTS: 

6 UNITS SCRATCH FILE 

7 UNITS, ANOTHER SCRATCH FILE OR TERMINAL 

THIS PROGRAM WAS DEVELOPED FRCM NOV. 1976 TO 

MARCH 1977 BY L. WOJTIS. 


XMAX=10. 

I=1 
INSEED=574657 
00Q(1) =1.0 
Q00(2)=2.5 
Q00(3)=25.0 
C0Q(4) =40.9 
000,(5) =50.0 
090(o)=60. 
000(7)=75.0 
Q00(8) =97.5 
Q00(Y) =100.9 


33 READ(5,30) NTI,IWP,KiK 


READ (5,431) TAREA, T3ASE, 1INC 
FORMAT (15,213) 
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GOs 4603. We Baily 319832) -KKK 
RZAD (5,833) Jour 
FORMAT (I 2) 

AIINC=IINC 

AREA=IAREA 

BMAX=0. 


vkd SYA. 


XNT=NTT 

DO 36 Id=1,200 
V (IJ) =00 

DO» 5,0. TLT=1),,50,0 
DO. 60 IZX=1, 20 
BIQITT,12X)=0. 
Adel IT, GX) R02 
FORMAT (15,5123) 
DO 51 ITT=1,20 
VA (ITT) =0. 

VB (ITT) =0. 
LPANTTs EDs 19 GO TO 31 


READ (8, 1, END=10) A(I,3) ,AS1,4S2,AS3 


DF s(An( Ly 3) wGb.2500.) G GO Oy 2 
FORMAT UES. 34/205) 
XMAX=AMAX1(X4AX,A (I,3)) 
I=I+1 
GO TO 2 
N=I-1 
vflsea Creat ( ZAP SMO) yas 
FORMAT, (8X, 4X,F7.2,A3,A2,A3) 
DO 701 I=1,N 
VIII) =A (1, 3) 
WRITE4/¢ 34) (Vid (1) ,T=1,8) 
IQI=0 
ZA=KMAX 
DO 702 I=1,N 


aF (VII (1) .LT.ZA) GO To 703 


GO TO 702 
TZE=1 

ZA=VII (I) 

CONTINUE 

IQI=IQI+1 
VdJJ {IQI) =VII (IZ) 

VII {IZE)=XMAX 
IF(IOI.GE. (N+1)} GO TO 705 
GO TO 704 


WRITE (7,37) (VJJ(1) ,I=1,5) 


DO 708 I=1,N 
Zitat 
VKK (I) =ZII*100./N 
WRITE (7,37) (VKK (1) -I=1,4) 
IVMI=N/JJIe 
FORMAT (1X, 2F 10.2, 3X,215) 


WRITE (7,710) VKK(IVM1) ,VKK{IVM2) ,LYM1, FVM2 


CALCULATION OF ALPHA AND BETA PARAMETERS FOR 


GAMHA DISTRIBUTION. 
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46 


SsUM1=0.0 
SUM2=0.0 
DO 150°T=1, 8 
SUM2=SUN2+AL0G (A (I, 3)) 
SUMN1=SUM1+A(T, 3) 
CONTINUE 
XM=SUM1/N 
XML=SUN2/N 
AA=ALOG (XS) -XML 
AX=1.04#4.0/3.0*AA 
AQ= (1.04#SQRT (AX) ) /4.0¥*AA 
BOQ=XN/AQ 
CQQ= BO*AD) 
WRITE(7,37) AQ,BO,COQ, XM 

SUM3=0. 

DO’ 732° Fn=aeu 


YUR=(A(IN,3) -XM) * (A (IN, 3) -XM) 


FORMAT (1X, 10F8.1) 
SUM3=SUM3+YUR 
YBN=N 
SUM4=SUM3/YBN 
SUMS=SUM4**.5 
WRITE(6, 3790) YBN 


WRITE (6,3790) SUM3,SUM4,SUM5 


SY=0. 

DO 4 T=1)K 
IRR=(A(I1,3)/5e) +16 

AI (IRR, 1)=AI (IRR, 1) +1. 
SY=SY+A(I,3) 

XN=N 

JX I= (XMAK/5~ 01. 

pO 500 S=tya x 

IXMP= (3-1) ¥*5 

Q=AT (J, 1) 100. /XN 

VYY (3) =¢ 

IZ=0*3. 

IF (1Z.EQs0)) GO TO 500 
CONTINUE 

SYM=SY/KN 

IW=2 

IWMP=IBASE 

WRITE (0,40) SYM 

IT=1 : 
TSYM=0. 

DO 18 IX=1,N 
XZZCLKpooe 

R(LX) =Oe 

SUMX=0. 

GO TO(46,47,48,609) , KKK 


KKK CHOOSES THE SAMPLING TECHNIQUE ACCORDING TO THE 
- ABOVE CODE. 


CALL GGTMAJ{INSEED,AQ,BQ, IWME,WWHK,K) 


FORMAT (1X, 2F8. 3,315) 


731 DO 14 I=1,1WMP 
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48 


609 


49 
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SUMX=SUMX+R (I) 
XZ2(1) =R EL) 
CONTINUE 
GO TO 4g 
ISTOP=LWNP/JJUP 
IVMN3=N-1VM1¥*JJIP 
XCIVMi=IV¥1 
£P=1,07x5CT 7" 
DOU? 13 -Texp=1,JIuP 
pO° 712 I=1,ESTor 
CALL GGU3(INSEED, 1,RX) 
TREEX/ AP 
IR=IR¢IVM1* (TEND- 1) 
XXIR=RX/AP 
XXIR=XXIR+IVM1¥* (IEND-1) 
CIR=IR 
YYIR=XXIR-CI2 
ICOQ= (JEND-1) *3+T 
PREY YR SLE .SORy GO To 729 
TTIR=YYIR/VKK (1) 
UDIR= (Vdd (I24+1V M341) -Vdd (TR+IVA3) ) *TTIR 
SUMX=SUMX+#VJ J(IR¢IVM3) tUUIR 
XZZ (ICQ) =VdU (IB+I V3) +UUIR 
GO TO 712 
SUMX=SUMX4+VJd (IR+IVMQ) 
XZZ(ICQ)=VdI (TR+IVES) 
CONTINUE 
CONTINUE 
GO TO 49 
DO 52 I=1,1WMP 
AP=1.0/XN 
GALL GGUS (INSEED, 1,8%) 
IR=RX/AP 
TF (R28) «GE.0.) GO TO S55 
R(IR)=1. 
SUMX=SUMN X4+A 4 TR, 3) 
KZZ (1) =4 (18,3) 
CONTINUE 
GO TO 49 
DO 53 T=1,7WMP 
AP=1.0/XN 
CALL GGU3(INSEED, 1,RX) 
UP (R LE Cs0) 1 o00 FO 4 
IR=RX/AP 
SUMX=SUMX+A (TR, 3) 
AZZ e Nene s) 
CONTINUE 
WRITE (6,37) (X22 (1) -1=1, LEMP) 
AMS=IWNP 
ASYM=SUMX/ANS 
W8 (IT) =ASYM 
BMAX=AMAX1(CMAX,ASYM) 
TSYMN=TSYMtASIY 
TSYM=(ASYM/5.) +1. 
CEUID) S25 Ye 
AL(ISYM, 38) =ALCISYM, IW) +1. 
TT2T1+1 
PF (LI. EE. NTI RCCOTION1? 
VD (18) =TSYM/XNT 
I8=IT-1 
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240 
241 
242 
243 
244 
245 
246 
247 
248 
249 
2590 
251 
22 
253 
254 
255 
256 
257 
258 
259 
260 
261 
262 
263 
264 
265 
266 
267 
268 
269 
2706 
271 
DZ 
273 
274 
AUS 
276 
207 
278 
279 
280 
281 
282 
283 
284 
285 
286 
287 
288 
289 
290 
291 
292 
293 
294 
295 
296 
297 
298 


650 


936 


934 


33:2 


931 


933 


935 


940 


99 


34 


102 


Wiles Gi, Se?) Bie Wi) 
VAR=0. 
DO 650 ICW=1,NTT 
VAR=VAR+ (K5 (ICK)-VE (Te) ) * (HE (ICW) -VB(I¥) ) 

VA (TW) = ((VARS/XNT) *¥*. 5) 

WRITE (6,237) VA(Ik) 
IYMP=IWMP+IINC 

ICH=1 

ICF=1 

QOPNT=NTT 
ICG=QOPNT*000 (ICH) J100. 
CA=XMNAXK 

DO 93h Sr 1 ANE 
IF(GEPEVSET SCA) GO TO 332 
GO TO 931 

ICE=I 
CA=CE(T) 
CONTINUE 

DE(TEPASOALGG) AGO TO, 933 
ICPF=ICF+1 

CE (ICE) =XMAX 
GO TO 934 
CF (ICH, IW) =CE(ICE) 
ICH=ICH+1 
IF (LENS GE. 10) GO £0. 93:5 
GO TO 9336 
IW=IW+1 

DO 940 T=1,5000 
CE(1I)=0. 

WRITE(7,300) IWMP 

IF (IWNP.GT.IWP) GO TO 99 
GO TO 100 

WRITE (7,34) NIT 

IZP=IW-1 
FORMAT(1X,15HNUMBRER OF DRAWS,18) 
JMP=(BMAX/5.) +1. 
IXPP=IwW-1 

FORMAT (1X,14,3X,20F6.0) 

DO 1058rJ=17,. 122 

BI(1,7J) =A1{1,1d) *100./XNT 
DO 104 IIl=2,JMP 

By (ete tai) = Gin ees) SOO / XN Ts ( Bat bt — ly liad) 0) 
CONTINUE 
CONTINUE 

FORMAT (15) 


PLOTTING THE PERCENTAGE FREQUENCY OF OCCURRENCE 
GRAPH FOR IMPACT ENERGY VALUES> 


CALLS PLOTN0L 0, Ue, 3) 
CALL PLOT (4.0,4.0,-3) 
CALL AXIS2(0.0,0.,'ZNERGY DENSITY IN JOULES/M*M',-28, 
«THO, OopOLPNO0. , 1) 
CALL AXIS2(0.,0.,'PERCENTAGE PREQUENCY OF GCCURRENCE' 
Ea Th 57'S) 0) ap Cho. 7 UNG G pee) 

CALI PLANS. ,3.5,5,NTT,AS1, AS2,A8S3,J3I5 2) 


TS 
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4AthD 
Ca -tp ALES I1K2 


Re 
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Oe 


_ 


= is 


Zar 
br : 


hd 


ee 


902 


901 


904 
993 


905 


9352 


PNR PIP Maya a) <e' 


114 


nO 920 I=1,438 

AX=1*25. 

ABC=AX/30 

CALL MDGAM(ABC, AC, PROB, IER) 
PP{T)=2ROB*100. 

PQ (IT) =AX 

PX (1) =PP(1) 

DOLIS 1 T=2,.48 

PX (1) =PP(I)-PP{I-1) 
IZA=1 

I=1 

SUE=0. 

pO 901 IJ=1,5 
IPP=I+IJ-1 

SUM=SUM+VYY (IPP) 
VRA(IZaA) =SUN 

T=It5 

IZA=IZA+1 

IF(IZA.GE.49) GO TO 904 
GO TO 902 

WRITE( 7,993) WRA(IK) ,IK=1,12) 
FORMAT (12F8. 2) 

CALL PLOL(O. 70 5,3) 

DO 905 I=1,28 

WHI=I 

WA=(WWI-1.)/4. 
VRB(I)=VRA(I)/10. 

CALL PLOT (WA,VRB(I) ,2) 
WB=WWI/4. 

CALL PLOT (WB,VRB(I),2) 
CONTINUE 

WA=.125 

U VAP), / 105 

CALL PLOT(WA,UV,3) 

CALL SYMBOL(HA,UV,.06,02,0.0,-1) 
DO 952 I=2,28 

WHI=I 

WA=(WWI/4.)-.2125 

De P Ketel 17, 10 

CALL PLOT (RA,UV,2) 

CALL SYMBOL(KA,UV ,.06,02,0.0,-1) 
CAGL PEGI = eA = iets) 


PLOTTING THE CUMULATIVE FREQUENCY GRAPH FOR 
IMPACT ENERGY 


CALL AXIS2(0.,0.,"ENERGY DENSITY JOULES/N*M',-25, 
+1 OR DE ne (Obes 1s) 
CALL AXIS2(0.,C.,'CUMULATIVE FREQUENCY! ,20, 
TORE Ue gone ges se ls) 
CALL PL1A(4.,3-5,5,NTT,AS1,452,483,JIIE) 
CELL BLOE IO. 0443) 
VC=VRB(1) 
DO 906 I=1,28 
WHI=I 
WA=HWI/4. 
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359 VM=VC/2. 


360 CALL PLOT (WA,V4, 2) 

361 CALL SYNBOL(WA,VN,.06,011,0.0,-1) 
362 VC=VCtVRB I+) 

363 906 CONTINUE 

364 GAM GPLOT (0. ,0.,3) 

365 DO 978 I=1,28 

366 WWI=I 

367 WA=WWI/4, 

368 UV=PP{I) /20. 

369 CALL PLOT (WA,UV,2) 

370 978 CALL »SYMBOL(WA,UV,.06,02,0.0,-1) 

371 € 

372 Cc 

3573 Cc 

“374 Cc PLOTTING THE STANDARD DEVIATION FOR A SAMPLE SIZE. 
3715 € 

376 Cc 

Sey (oe 

378 37 FORMAT (1X, 10F8. 2) 

SYS) CAO) ICOR (Aa As pas) 

380 GALE VEROGAX (0. 70, NUMBER! 7-6, 

381 Hsia Qiere 0s 9-20 20's) 

382 CALL LOGAX(90.,0.,'"STANDARD DEVIATION', 18, 
383 6 7 AO KO eZee) 

384 CALE MEE WAs( Dies, De gp hKK p ND Lg AS le ASS ¢ ASS, OU £) 
385 ARA=IBASE 

386 I=2 

387 908 EPO MRA! 20.) ORs (YA (1) <LE..0.).) GO TO 9381 
388 AVT (I) =ALOG10 (ARA) 

389 AVB(1) =ALOG10(VA({I)) 

390 98% ARA=AKAtAIINC 

391 WRITE (6,37)ARA,VA(I) ,AVT(I) ,AVB UTI) 
392 I=I+1 

393 TE tLeGT.£72) GO .PO,907 

394 GO TO 908 

39/5 907 NOP=IZP 

396° & CAELAPEOR (IO. 0.23) 

397 AQP=AVT (2) *5./2. 

398 AQR=AVB(2) *5./2. 

399 GALL ESYMBOL( AOR, NOR, .06, 911, 050),—1) 
400 CALL SHOOTH (AQP,AQR,0) 

401 IXZ=NOP-1 

4Q2 DO 909 I=3,1XZ 

403 AQP=AVT (I) *5./2. 

404 AQR=AVB (I) *5./2. 

405 CALL SHOUTH (AQP,AQR,-2) 

406 909 CONTINUE 

407 AMM=AVT (NOP) *5./2. 

408 AMN=AVB (NOP) *5./2.6 

409 GALE SHOOTH (AMM, ANN, -25) 

410 DO 955 I=3,NOP 

411 AQP=AVT({I) *5./2. 

412 AQR=AVB(I) *5./2. 

413 955 CALL SYMBOL (AQP,AQR, -06,011,0.9,-1) 
414 GALL SPLOT (0..,- 15. ,-3) 

415 Cc / 

416 C 

417 Cc 

418 c PLOT OF THE IMPACT ENERGY POR A GIVEN SAMPLE SIZ” 
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444 
442 
44 3 
444 
445 
446 
447 
448g 
449 
450 
451 
452 
453 
454 
455 
456 
457 
458 
459 
460 
461 
462 
463 
464 
465 
466 
u67 
468 
469 
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471 
472 
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834 
835 


a2 
911 


JH7/ 


956 


836 


839 
838 


841 


CALERUOGAMHO. 0. , (NUMBER ,-6,5.0,0., 00 ,-2e,0-) 
IF(KKK.EQ.4) GO TO 834 
CRig, PAT SA(0.g 6. CENERGY DENSITY® 14,10. ,90.,.00, 20 0,-1.2) 
GO TO 835 
CAMs AXUS21006> 0.7, "ENERGY DENSITY", 14,10. ,90.¢0:,10.,-1.) 
GALIHSPEOT (05,0. , 3) 
CRE) PRIA (To, 4 «, KKK, NTT, AST, AS2,AS3,JJ0P) 
IP(KKK.EQ.4) GO TO 836 
DO 911 I=2,38 
AXX=IBASE 
CQ=CF( E,:2)'/20, 
BXX=(ALOG10 (AXX)) *2.5 
CALLSSYUSOL (EXX, CO,2 06, 1,0.0,-1) 
DO 912 J=3,1ZP 
AXX=AXX+ATINC 
CQ=CF (I,J) /20. 
BX X= (ALUG10 (AXX)) *2.5 
CALLE SYMBOL CBX XACO. 06, li) «Og = 1) 
CONTINUE 
DO“ 956 i=2;8 
AXX=IBASE 
CQ=CF (1,2) /20. 
BXX= (ALOG10 (AXX)) *2.5 
CALL, SMOOTH, CEXXU.CO,,0) 
DO 957 J=3,1XZ 
AXX=AXX+AIINC 
CQ=CF (I,J) /2C. 
BXX=(ALOG10 (AXX)) *2.5 
CALL SMOOTH |BXX,CQ,=2) 
AXX=AXX+ATIINC 
CQ=CF(I,NOP) /20. 
BXX=(ALOG10 (AXXK)) *2.5 
CALLe-ShOOT Hy bX y).CO, 125) 
CONTINUE 
GO TO 837 
DO 838 T=2,8 
AXX=IBASE 
CORCPADIZ} / 1.0. 
BX X= (ALOG 10 (AKX) ) *2.5 
GALTS SINBOLN( EX, Col, 06, 2,0. 0,—1) 
DO 839 J=3,12P 
AXX=AXX+AIINC 
CQ=CF(I,J)/10. 
BYX=(ALOG10 (AXX)) *2.5 
GALINSVEDOIe BEX, CO, .06,15,0.0,-1) 
CONTINUE 
DO 840 I=2,8 
AXX=IBASE 
CQ=CF(I,2)/10. 
BXX=(ALOG10 (AXX)) *2.5 
CALL SMOOTH (BXX,CQ,0) 
DO 841 J=3,1XZ 
AXX=AXX+AIINC 
COQ=CF(I,J)/10. 
BXX=GALOG1C (AXX)) *2.5 
CALL SMOOTH (Bi KOO, 2) 
AXX=AXX+ALINC 
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498 
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504 
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506 
507 
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509 
510 
511 
St2 
Sts 
514 
S1S 
516 
S47 
518 
519 
520 
521 
S22 
523 
524 
525 
526 
527 
528 
529 
530 
531 
532 
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535 
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537 
5338 


840 
837 


991 


ALOLOLOLO Orne 


9153 


954 


359 
958 


CO=C 
BX X= ( 
CALL 
CON 
CALL 
Do 99 
UG=T 
UF=7. 
URR=U 
UH=00 
CALL 
CALL 
CALL 
CALL 
CALL 


PLOT OF THE PROBABILITY OF OCCURENCE FOR AN IMPACT 


ENERG 


GALL BAX S260 Ue, "ENERGY IN JOULES/M*M* , 


F220 7 


F(I,NOP)/10. 

ALOG10 (AXX)) *2.5 
SMOOTH (BXX,CQ,-25) 

SNE 

PLODY Bie igi7:e ty 3) 

1 I1=2,8 


$0—(UG=1) *.125 
F-.1 

Q (1) 
NUMBER(8.1,US,.06,0H,0.,1) 
SYMBOL (82 5,UP,.06,1,0-,—-1) 
SYMBOL(8.6,UF,.06,1, 0«,-1) 
SYMBOL(8.7,UF,~.06,1,0.,-1) 
PLOMU1S%. » 0%..73) 


Y VALUE. 


op Gop Wan2 Se py Ue) 


GALL WAGES 2(05,0., {PERCENTAGE FREQUENCY (OF OCCURRENCE”, 


*34,7. 


CALLE, yRL WA (Get so pkKKK NDT ,AS1T,ZASZ,AS3,JIJE) 


DO 95 
CALL 

XCP=A 
CALL 

DO 95 
XOP=I 
XOP= ( 
XCP=A 
CALL 

CONTI 
XCP=A 
CALL 

CONTI 
DOLY5 
CALL 

XCP=A 
CALL 

DO 95 
XOP=I 
XOP= ( 
XCP=A 


CALL S¥MBOL(KOP,XCP) .06,0,0.,7-1) 


CONTI 
CONTI 
CALL 

UH=IB 
DO 96 
UG=I 

UF=7. 
URR=U 


CALLYSYEBOL (Ss 9, UF ,.06, 30K =,0.,3) 


CALL 


POA Oo 5 OK Aas.) 


4 4J=2,12P 
AECURMOS AUS Ash! 
P(Mp Ie 10, SANT 
SMOOTH (.1,XCP,0) 
3 1F=2,, 34 


XOP*.2)-.1 

I (1,5) *10./XNT 
SMOOTH (XOP,XCP,-2) 
NUE 

1(35,d) *10./KNT 
SMOOTH (6.9,XCP,-25) 
NUE 

avd=2,E2P 

PEOT (09,09, 3) 

Ld gd)* 105 /XND? 
SYUBO(.1,X%CP,. 06,5, 06--1) 
Jota? ,35 


XOP*.2)-.1 
I(I,J) *10./XNT 


NUS 
NUE 
PLOT (8.,7.,3) 
ASS 


OM E=2 1-22 


i321 UG> 3.) *..125 
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960 
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CALL “SYADOL (6. 5,UF +. 06,1, 005-1) 
CALL SYMBOL (6.6, UF,.06,1,0.,-1) 
UH=UB+AIINC 

CALC SYPCOLL0. L.U Fn. 0651¢0., 1) 
CALL PLOT(0.,15.,-3) 


PLOT OF THE PROBABILITY OF OCCURENCE FOR AN IMPACT 


ENERGY VALUE. 


CAV WAX S20 cbs SONERGY EN JOULES /M*M® , 
SVL Oi ieee AON AO SKe. vileow) 


CALL AXIS2(0.,0.,'PERCENTAGE FREQUENCY OF OCCURRENCE', 


SU ee OG HOw Oe ay 


DO 988 S=2,IZP 
CALL, PL TAOO sus Oleg 3) 
CALL, PL AA Gene do DnkKK NTT, AS 1 AS2,4S3, 3002) 
XCP=AI(1,J) *10./XiT 
CALL SYNBOL({.2,XCP,-06,3,0.,-1) 
DO 989 T=2,35 
XOP=I 
XOP=XOP*,2 
XCPHAT (1, J) ¥ 10. /XNT 
CALL, SYMBOL (XOP., SCP pus 6504.0 -,—-1) 
CONTINUE 
CONTINUE 
CALL SLOT (Oe « sa era) 
UH=I BASE 
DO 380 I=2,12ZP 
UG=I 
UP=7.13-(UG-3.) *.125 
URR=UF-.1 
CALL SYMBOL (5.9,UF,.06,3HK. =,0.,3) 
CALL NUMBER(6.1,UF,-06,UH,0.,1) 
CALL SYABOL loro, Ur, 5 06,1,0.,=1) 
CALL SYU50Ldb. 6, UF ys 06,1;0.,-%) 
UH=UHtAIINC 
CALL SYMROL(6.7,UF,-U6,I1,0+,-1) 
CRUEL PLOT 412.5 0a <3) 
CALL LOGAX(0.,C.,*NUKBER PER 1000 KM¥KM!,-21, 
#58 ei eee, abla! 
CEELMPLIA (Gs 1425 7RRK, NTT, AST, AS2,AS3,IIIP) 
CALL LOGAX(0.,0.,'STANDARD DEVIATION', 18, 
es ee ee oa Oa) 
XAREA=1000. /ABREA 
ARA=IPASE 
T=? 
AVT (1) =ALOG10 (ARA*XAREA) 
ARA=ARA+AIINC 


I=I+i 

Tee te i). Gh) Te) yao? 
GOTO .858 

NOP=IzZ2 


CREP LOT Cs 0s. 3) 
AQP=AVT (2) *5./2. 

AQR= (AVZ(2) +1.) *5-/2- 

CALL SYMBOL (AQP, AOR, -06,011,0.0,-1) 
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599 CALL SMOOTH (AQP,AQR, 0) 


600 IX Z=NOP-1 

601 DORS OIF r=Si Ke, 

602 AQP=AVT (1) *5./2. 

603 AQR= (AVB{1) +1.) *5./2. 

604 CALL SMOOTH (AQP, AQR, -2) 

605 809 CONTINUE 

606 ANM=AVT (NOP) *5./2. 

607 AMN= (AVB (NOP) +1.) *5./2. 

608 CALL S“OOTH (AMM, AMN,-25) 

609 DO 855 I=3,NOP 

610 AQP=AVT (I) *5./2. 

611 AQR= (AVB (1) #1.) *5./72. 

612 855 CALL SYNBOL( AQP, AQR,.060,011,0.0,-1) 
613 CAED PEOT (CG. ,- 15. 5-3) 

614 CALLVCEOGAX (0076.6 7 UNU MBER PER 110100 *KRMSKM 2 7=21, 
615 3 AWls GUcR Sea AWa)) 

616 CALLS LOGAX (0170.6 7 UME AN TDIBEE RENGE BR ERWEEN CEST EMA TES % 7 
617 EO ll G57 NOG pS liege earls) 

618 CALLEPL1Ng6s 24 55, KKK, NIT, AST AAS] ZAS3 7IIIB) 
619 XAREA=1000./AREA 

620 ARA=IBASE 

621 I=2 

622 256 AVT{1I) =ALOG10(ARA*XAREA) 

623 WRITE(6,37)AVT(I) ,ARA,VA(I) 

624 ARA=ARAtATINC 

625 I=I+1 

626 TF (lL. GTeiZP)) GO TO 257 

627 GO TO 256 

628 2S) NOP=I2P 

629 CALLER LOM iOeng Olenso) 

630 WRITE(6,37) XN, (VA(IVMP) , IVMP=2, NOP) 
631 AQP=AVT (2) *2.5 

632 AQR=(VA(2)) *100./VB (2) 

633 261 COR= (ALOG1C (AQK) +1.) *#2.5 

634 WRITE(6,37) AQR,COR 

635 GALLS YMBOMGAOP, GOR 160.6% Olt, 0s10)¢ — 1.) 
636 CALL SMOOTA(AQP,COR, 0) 

637 IXZ=NOP-1 

638 DO 258 J=3,1XZ 

639 AQP=AVT (I) *2.5 

640 AQR= (VA(I)) *100./VB (1) 

641 263 COR={ALOG10 (AQR) +1.) *2.5 

642 WRITE (6,37) AQR,COR 

643 CALL SMOOTH {AQP,CQR,-2) 

644 258 CONTINUE 

645 AMMN=AVT (NOP) *2.5 

646 AQR= (VA (NOP) ) *100./VB (NOP) 

647 265 COR=(ALOG10 (AOK) +1.) ¥2.5 

648 WRITE (6,37) AQR,CQR 

649 CALL SVOOL NN CAM; COR, —2 ©) 

650 DOM 27 Mae), kOe 

651 AQP=AVT (I) *2.5 

652 AQR=(VA(I)) *100./VB (1) 

653 273 CORFU OGMOIA OR) +1.) = 2.60 

654 PRY CRELMSYUBOLIINO@L COR, © 06, O11, O00 5—1) 
655 CALL PLOT{10.,0.,999) 

656 SUOye 

657 END 


658 SUBROUTINE PLIA(XST,YST,KKK,NTT,AS1,4S2,4S3,J3JJP) 
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659 
660 
661 
662 
663 
664 
655 
606 
667 
668 
669 
670 
671 
byhZ 
673 
674 
675 
676 
677 
678 
679 
680 
681 
682 
683 
684 
685 
686 
END OF FILE 
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120 


SUBROUTINE TO PLACE THE DATE OF STORM ON THE VARIOUS 
PLOTS. 


CALL’ SYMBOL (XST, YST, . 107A51,0.,3) 
CALL™SYHBOL (XST+.4,YST,.10,AS2, 0. ,2) 
CALL SYMBOL(XST+.6,YST,.10,A453,0.,3) 
UST=YST-.4 
GO TO (251,252, 253,254,256), KKK 
CALL SYMBOL(XST,UST,.10,19HFITTED DISTRIBUTION,0.,19 
GO To 255 
CALL SYMBOL(XST,UST,.10,19HWITHOUT REPLACEMENT,0.,19) 
GO TO" 255 
CALL SYMBOL(XST,UST,.10,16HWITH REPLACEMENT, 0,16) 
GO TO 255 
CALL SYMBOL(XST,UST,.10,19HSTRATIFIED SANPLING,O.,19) 
CALL S¥rp0tusss, UST-.6,.10,2HM=,0.,2) 
XPJJ=JJU? 
CALL NUMBER(XST+.5,UST-.8,.10,XPJJ,0.,-1) 
CALENSYEBOR(XST, UST-.4,.10,2HN=,0.,2) 
XNIT=NTT 
CALL NUMBER (XST+.5,UST-.4,.1C, XNTT,0.,-1) 
RETURN 
END 
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NEEL TCAn hOmeck MNONs  GAr LO) Us hep TO 
HAILPAD IMPACT ENERGY VALUES 
DIMENSION A(17120,3),8 (200) ,AI (120,2C) ,V(250) 
DIMENSION VA(20),B1(120,20) ,VB (20) ,CF (9,20) 
DIMENSION VYY (50C),VX (60) ,VY (60) ,CE(500C) 
DIMENSION AVT(2C) ,AVB(29) ,YRA(60) ,VRB(60) 
DIMENSION PX (60) ,PQ(60),PP(60) ,OOQ0(Y) 
DIMENSTON PPP (120) ,2OP (1290) 7 IPMS (2) 
DIMENSION VII(120) ,VJJ (120) ,VKK (120) 
DIMENSTFON TV XRXK (SOY SV YYYY (50) , VIT (10) 
COMMON #W(30,30) ,W8(200G) 
COMMON BUF (2048) 
CALL PLOTS (BUF,8192) 
XMAX=10. 
I=1 
INSEED=123457 
READ(5,30) NTI, I1WP,KKK 
READ(5,431) IAREA 
READ (5,432) IZK 
FORMAT (12) 
FORMAT (15,213) 
AREA=IAREA 
BMAX=0. 
XISY=0. 
XNT=NTT 
I=1 ; 
READ (8,1, END=C4)A (1,3) ,AS1,AS2,AS3 
FORMAL (CS XOX Es 3g AZ, AS) 
I=I+1 
GO TO 3 
SY=0. 
N=I-1 
DO 23 I=1,N 
SY=SY+A (I, 3) 
XN=N 
SYM=SY/XN 
DOSS 0Mdr l= 50'0 
DU 50 IZX=1,20 
BI (ITT, 1ZX) =0. 
AL (ITT, 1ZX) =0. 
ICQ=1 
DMS=0.0 
pO 611 I=1,30 
DO 611 J=1,30 
WW (I,J) =0.0 
READ(1;, 612, END=613)18,19,AI8 
FORMAT (22 AX, 12, 7% F152 
WW(I8,19) =AI8 
GO TO 614 
PO RMA (eu or, 8s) 
DO Sy) RAE AAO, 
VA (ITT) =0. 
VB (ITT) =0. 
FORMATCVX, 2F 1022,.3X%,215) 
TwW=2 
IT=1 
VAR=0. 
TSYM=0. 
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DO. 16. 1%=1;'N 
Rq( IX) =0. 
SUMX=0. 
AMS=0. 
DO 644% IDQ=1,1C0 
CALL GGU3(INSEED,2,R) 
XSUM=0. 
XNN=0. 
I1=1 
J1=1 
TFAGUUH GL wil). bO.O.) OR. (WW(T1,01+1ZK) .EO0.0.0)) GO T: 
Trower gLytlZ*,J1).6O.C.)-.OR. (WH(IL1I+IZK,J1+IZK).EQ.0.)) Gt 
WAI=WH (11,01) -NW(11,0141ZK) 
LF (WA 1 a BevO0,) 7 GHA T=— 164A 
Vi=WW (11,01) +WA 1*R(1) 
WA2Z=WW(I1+1ZK,J1)-ew(11+IZk,J1+1IZK) 
IF(HA2.LE.0.0) WA2=-1.¥*WA2 
V2=WW(IT1+1ZK,J51) +W¥A2*R (1) 
VM=(V1+V2)/2. 
WA3=WH(I1,31)-NW(I1+1ZK,J1) 
IF(WA3.LE.0.0) WA3=-1.*9A3 
W1I=RW (11,01) +#A3*R (2) 
WAG=WW({I1,014+12K)-WH(T1+IZK,J1+1IZk) 
IF(WA4Y.LE.0.0) WA4Y=-1.¥*W AY 
W2=WW(L1,01+1IZK) +8A4*R (2) 
WM=(H1+W2) /2. 
VAL=(Vitit+WM) /2. 
XSUM=XSUM+VAL 
XNN=XNN+1. 
FORMAT Hix 4 2ees;, JE Ye 2) 
J1=J 1+IZK 
EFAS IsiGh229) GO TO 5 
GO TO 7 
J1=1 
Ti=1I1+1ZK 
IFgdils GES 29) 4GO/TO 8 
GO TO 7 
SUM X= XSUM+SUHX 
AMUS=XKNN+AMS 
ASYM=SUMX/AMS 
W8(1IT) =ASY 
QCI=ICQ 
IWMP=XNN/OCI 
TINC=XNN/OCI 
TSYS=TSY4+ASYM 
TSYHSMASY 50) #0, 
CEC(IT) =ASYM 
AI{ISYM, IW) =AL(ISYM, 1W) +1. 
IT=1T+1 
ITF (AT. LEGNTT) 4 COPTO 17 
VTT (ICQ) =AMS 
VBq1I¥) =(TSYA/XNT) 
I8=IT-1 
WRITE (6,37) VS(19) -AMS, XNN 
YVAR=0. 
DO; 650 pF CUH4, NTT 
VAR=VAR#t (88 (I1CW) -VB (IW) ) * (#8 (TCW) -VB (I) ) 
VA(IW)=( (VAR/XNT) **.5) 
WRITE(6,37) VA(1Zk) 
ICQ=ICQ+1 


‘T Co Meo Os O% . (STAIE, 
Fs) {tere Ga. CHST# 10 8 


4 


‘ rt ly (Ean STE 


get rie 
a¥~ (WOT) AW) +AAV=MAY 
J=HCULPAY 
Ty AV (VEO) STIRa 

Me ital 


(AN) oY (93) 849 i 


= 
> 


S e-enes 


oa, 


_ 


. f - 
et , s = 
; aD =. Wes 5 = _ 
7 i? q ‘ke % i - 
- - - a. . = a ee 
. : z - | ihe 
oe. is » ‘» = - 
4 7 s: ‘ 
‘ s—. : 
- a 3 = ‘ . - = - 
= a ow z _ a ave - 7 
. = ¥ —. os er <2 
i - - > = a ” x 
P we q ‘ " 2 4 ~ : 
a - a me be , : / = 
= oo = f al = ho ‘y ¢ 
7 * ~ > , Z 


j 
, é 
af, 
=f 
, D 
thie 
¥ 


335 
SS i/ 


99 


34 


102 
300 
31 


908 


907 


909 


955 


QR PIAA? QD IM) 


IW=IW41 
FORMAT (1X, 10P8. 2) 

WRITE (7,300) IWMP 

IF(ICQ.GE.4) GO TO 9¥Y 

GO TO 100 

WRITE(7,34) NTT 

IZP=1i-1 

FORMAT (1X, 15HNUMBER CF DRAWS,T8) 
JMP=(BMAX/5.) +1. 

IXPP=IW-1 

FORMAT (1X,14,3X,20F6.0) 

FORMAT (15) 

CRLEUPLOT (0.050% ,3) 

GALE / ®L074450,4.50,-3) 

CALL LOGAX(0.,C.,'NUMBER',-6, 
£505 O0n71.,-2. 0). 

CALL LOGAX(0.,0.,"STANDARD DEVIATION', 18, 
Ore ¢ Oa Oe sce Os) 
CALE (Pid An sou bo pKKK, NTT, AS 1, AS2,AS3) 

JPP=1 

ARA=VTY (JPP) 

T=2 

AVT(I) =ALOG10 (ARA) 

AVB(1I) =ALOG10 (VA (I) ) 

WRITE (6,37) ARA,VA(I),AVT(I) ,AVB(T) 
I=I+1 

JPP=JPP+1 

ARA=VTT(JPP) 

IF(I.GT.IZP) GO To 907 

GO TO 908 

NOP=IZP 

GALL PLOT (Oe 0s 72) 

AQP=AVT (2) *5./2. 

AQR=AVB(2) *5./2. 

CALL SYMBOL(AOP,AQR,-060,011,0.0,-1) 
CALL SAOOTH (AQP,AQR, 0) 

IXZ=NO?-1 

PO 9099t=3 7 XZ 

AQP=AVT(I) *5./2. 

AOR=AVB(I) *5./2. 

CALL SMOOTH (AQP, AQK,-2) 

CONTINUE 

AMM=AVT (NOP) *5./2.6 
AMN=AVB (NOP) *5./26 

CALL SHOOTH(AMM,ANN, -25) 

DO 955 I=3,NuP 

AQP=AVT (I) *5./2. 

AQR=AVB(I) *5./2. 

CALL SYMBOL(AQE,AQR, -06,011,0.0,-1) 
EALEP PLOT (ON, 15 e473) 


PLOTROLetHh SRORBALLLITY OF OCCURENCE FOR AN EM PACT 
ENERGY VALUE. 


CALL AXIS2(0.,0.,'ENERGY IN JOULES/M¥*M', 
Bea O fe plicig Vie ena iss) 
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CALL AXTS240.,0.,'PERCENTAGE FREQUENCY OF OCCURKRENCH", 
* 34,1. SONOS ale scie) 


CALL. PLIDA( Geip4 > So KKKG NTE, AS1,AS2, 453) 
DO. 954. Jz2,12P 
CELL PLOT (02,0. 93) 
XCP=AI (1,J) *10./XNT 
CALL SNOOTH(.1,XCP,0) 
DO 953 I=2,34 
XOP=I 
XOP=(XOP*.2)-.1 
XCPSAI 41, J) 340. /XNT 
CALL SMOOTH (XOP,XCP, -2) 
CONTINUE 
XCP=AI(35,d) €10./XNT 
CALL SMOOTH(6.9,XCP,-25) 
CONTINUE 
DO 958 J=2,12P 
CALL PLOT (0. .0. #3) 
FCPSRI EI .d)% 70. ARNT 
CALL SYMBOL(.1,XCP,.06,d,0.,-1) 
DO 959 I=2,35 
XOP=I 
XOP=XOP*, 2 
XCP=AI (I,J) *10./XNT 
CALL SYMBCL(XOP,XCP,.06,J,0.,-1) 
CONTINUE 
CONTINUE 
CHILL LPLOT (8. 57573) 
JPP=1 
DO 960 I=2,I2P 
UH=VTT {JPP) 


UG=I 
UR=7... 1d stUiG=s) *. 125 
URR=UF-. 1 


CALL NUMBER (5.9,UF,.06,3HK =,0.,3) 
CALL NUMBER (6.1,UPF,.06,UH,0.,1) 
CALL SYMBOL(6.5,UF,-06,1,0.,-1) 
CALL SYMBOL(6.6,UF,.06,1,0+,-1) 
JPP=JPP+1 

CRE. cS VAL > Jreill Fee Oils Oe = 1) 
CALL PLOT (.15.,0.,-3) 


PLOT OF) THE. PROBABLLITY, OF OCCURENCE FOR AN IMPACT 


ENERGY VALUE. 


CALL GAXGUSZID eyes VENERGY EN JOULES /ItN, 


OVX OT Are 0) aust APSE pao) 


CALLA AKESZ (On, Ue RE ERECENTAGS FREQUENCY OF OCCURRENCE, 


«30 ee Wi | Oan ~ Ae) 
DO, 98p J=25 lee 
CALE OCOD S fis) 
CALL PL athe de ab, KKK, NTT,AS1T,AS2,AS3)} 
XCP=AI (1,J) ¥19./XNT 
CALL (SY Sa(s2¢kCP as 06,0, 0.571) 
DO 989 I=2,35 
XOP=I 
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XOP=XOP*. 2 
XCP=AI (I,J) *10./XNT 
CALL SY UBL (xOR CE, 06,0900 »- 1) 
CONTINUE 
CONTINUE 
ExT SeLuUl Ce ieiesh 
JPP=1 
DO 980 I=2,IZP 
UH=VTT (JPP) 


UG=I 
UP=7.13-(UG-3.) *.125 
URR=UF-.1 


CALL NUMBER (5.9,UF,.06,3HK =,0.,1) 

CALL NUMBER (6.1,UF,.00,UH,0.,1) 
CALL SYUBOL (6..5,,0F.. 06, L, 0. ¢=4) 
CALL SYMBOL(6.6,UF,.06,I,0.,-1) 

JPP=JPP+1 
CALL SY4BOL(6.7,UF,.06,1,0.,-1) 
CALL PLOT (0. ,-15. ,-3) 

CALL LOGAX(0.,0.,*NUNBER PER 1000 KM*KM',-21, 
SOR BOK plays Bees) 

CALE (SE 1N6 uy 4S KKK, NTC,AST,AS2,NS3) 
CALL LOGAX(0.,0-,*STANDARD DEVIATION', 18, 
AA ASG HOU a pees OG) 

XAREA=1000./AREA 


JPP=1 
ARA=VTT (JDP) 
I=2 
AVT (I) =ALOG10 (ARA*XAREA) 
Tai) 


FF(1.Gloiz Py GOuTO. 657 

JPS=IPP+1 

ARA=VTT (JPP) 

GO TO 858 

NOP=IZP 

CRLEL, PLOT(0.,0.,3) 

AOPEAVT (2) *5./2. 

AQR= (AVB(2) +1.) *5./2. 

CALL SYMBOL(ACP,AQR,-06,011,0.0,-1) 
CALL SMOOTH {AQP, AOR, 0) 

IXZ=NOP-1 

DO 809 I=3,1XZ 

AQP=AVT(I) *5./2. . 

AOR= (AVB(1) +1.) *5./2. 

GALL SMOOLH (AO2, AOR, —2) 

CONTINUE 

AMM=AVT (NOP) *5./2. 
AMN=(AVB(NOP)+1.) *5./2. 

CALL SMOOTH (AMM, AMN,~-25) 

DO 855 I=3,NOP 

AQP=AVT (I) *5./2. 

AOR= (AVB(1I) +1.) *5./2- 

CALL SYMBOL(AOP,AOR,-06,011,0.0,-1) 
CREEP ROTC) 527.0 6.53) 

CLT TOGAME(Ome, On UNO US ER BIER 1HOGC KMS Kiet = 24, 
SO Als pO) Sh A Sed AO 

CALEMLOGAXK(Ue Oe USAN DEFFERENCE  BETWERN ESTIMATES * 
OS Uo Dy M05 p= done sian OS) 
CALL PEMACO. 5.4. OD, KKK, NTT, ASL, ASZ,AS3) 
XAREA=1000./AREA 
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JPP=1 
ARA=VT? (JPP) 
I=2 
AVT (I) =ALOG 10 (ALA*XAEREA) 
I=I+1 
IF (MAGISIWS) GO To 257 
JPP=JPP+1 
ARA=VTT (JPP) 
GO TO 256 
NOP=IZP 
CRUE CPLOT (0% 0007 3) 
XM=SYM 
WRITE (6,37) XM, (VA(IVMP) , IV4P=2, NOP) 
AQP=AVT (2) *2.5 
AQR= (VA (2) ) *100./V23 (2) 
CQR=(ALOG10 (AQR) +1.) *2.5 
WRITE (6,37) AQR, COR 
GALL SYMBOL CAO2, GOR, .06,011,.0. 0,—1) 
CALL SMOOTH (AQP,COR,0) 
IXZ=NOP-1 
DOP 258) t=3), UNZ 
AQP=AVT(I) *2.5 
AQR= (VA (I) ) *100./VB (I) 
CQR={ALOGI10 (AOR) +1.) *2.5 
WRITE(6,37) AQE,COR 
CALL SMOOTH (AOR, COR, =Z2) 
CONTINUE 
AMM=AVT (NOP) *2.5 
AQR= (VA (NOP) ) *100./YB(NOP) 
CQR= (ALOG10 (AQR) +1.) *2.5 
WRITE (6,37) AQR,COR 
CALL SHOOTH (AMM,CQR,-25) 
DO 271 2=3, NOP 
AQP=AV?{I) *2.5 
AQR=(VA(I)) *100./VB (I) 
COR= {ALOG 10 (AQR) +1.) *2.5 
GAMING FS orl SO (MOR COR a Oly Oleg Ore Olpi—) 1) 
CIMEIh ILO CUM oA Ula pO 
STOP 
END 
SUBROUTINE PLIA(XST, YST, KKK, NTT,AS1,AS2,AS3) 


SUBROUTINE TO PLACE THE DATE OF STORM ON THE VARIOUS 
PLOTS. 


CALL OSY REO UKs, VST, + 10, WS7, 0673) 
CALL SPapOlmsi+.4,15T,.1.0,852,0+92) 
CALEVS TMBOL WKS t eG, VST, . 10,8535 0%, 3) 
UST=YST-.4 
IF(KKK.EQ.9)GO TO 256 
GOTO. (2515297, 253,254, 256), KKK 
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GALES NIMBOLIGKST ,UST, «MO, 1 I Me lu beD DISTREBUT LON, 0. 512) 


GOTO; 255 


CALL SYMBOL(XST,UST,.10,19HWITHUUT REPLACEMENT, 0O., 19) 


GO LOZ 55 


CALL NOM OlCMol UST, . NO, Ou GLH REPLACE MEN 7 0l, 16) 
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359 GO) 202255 

360 254 CALL SYMBOL(XST,UST,.10,19HSTKATIFIED SAMPLING,O.-1%) 
361 (0) lO). BBS 

362 - 256 CALL SYMBOL(XST,UST,.10,11H2-D CONTOUR,S., 11) 

363 255 CALL SYSBOL (XST, UST--4,-10,2HN=,0.,2) 

304 XNTI=NTT 

365 CALL NUMBER(XST+.5,UST~.4,.10,XNTT,0.,~-1) 

366 RETURN 

367 END 


END OF FILE 
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DIMENSION AA (22,25) ,X (6CO) ,¥ (600) ,Z (600) ,CVAL(10) 


INTEGER IOP (10) /10*0/ 
COMMON A (8) ,BUF (2048) ,B(25) ,VOP(10) 
DATA AO/*YES*/ 
CALL PLOTS (SUP,€192) 
LOP (1) =1 
Top (2) =1 
Lop (3) =0 
IOP (4) =1 
NC=10 
CKIN=25.0 
CSET=CMIN 
CINT=50.0 
FORMAT (1X, 13P8.2) 
YEAR=1976 
XMAX=-5.0 
XMIN=15 20 
YMNAX=0. 
YMIN=22. 
IEN=1 
Tf 
READ (8, 1,END=10) (A (II) ,II=2,4) ,AI3 
PORMAT (2X,3F2.0,4X,7F 702) 
CALL GENER (JR) 
¥T=4. (2) 
XX=A (1) 
XMAX=AMAX1(XMAX,XX) 
XMIN=AMIN1(XMIN,XX) 
YMAX=AKAX1 (YMAX YY) 
YMIN=AMIN1V(YMIN,YY) 
X (IMM) =XX 
Y (IMM) =YY 
Z (INN) =AI3 
INM=INM+1 
I=1+1 
GO, TO 2 
N=1NM 
XSIZE=KMAX-XMIN 
YSIZE=YNAX-YMIN 
WRITE (6,310) XSIZE,YSIZE 
FOKMAT (1X,2P10.2) 
BO 2 0,0teal, BC 
CVAL (I) =CSET 
CSET=CSET+CINT 
CONTINUE 
CALL PLOT (0.,0.,3) 
DO 300 I=1,N 
VX=X (1) 
VY=yY (I) 
CALL SYMBOL(VX,VY,-06,3,0.,-1) 
CALL NUMBER (VX+.09,VY,.05,Z(I) ,0-,1) 
CONTINUE 


Chi eoGCONT(T Sl Ze. YS! Dols eels beGVhles 
aNC,IOP, VOP,XMIN,XMAX, (MIN, YMAX) 


CALL PLOT(0.,0.,999) 
STOP 
EKD 

SUBROUTINE GENER (KTR) 
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COMMON A(8) ,22(9) ,BUP (2048) 
A(1)=1. 
A (2) =21. 

IF(A(4) -E0.29) A(2)=23. 
N=(A (2) -1) fo. 

N=(A(1)-1) /2- 
YM=N 
YN=N 
R=A (4) 
T=A (3) 
Y=(6.* (T-19.) +YM+YN/2. +. 25) *1.6084 
IS=A (2) 
K=r4+1 
GO SOUL(27 (28/2728 727), 28), ,.K 
P=IS-— (14+6* (K-1) ) - 
GO TO 29 
P=K*6-IS 
IP(A(1)-L=.2) GO TO 30 
N=(A(1) -2.)/2- 
GO TO 31 
H=A (1) 72 
XN=N 
L=(A(3)-15.)/4. 
Z=1.3 
C=5i2 

Te (Le LEeeeye cour. 32 

Sie. 


TFP(R-LE.10)) GO TO 33 
X=(6.* (29.-R) #C-Z* (L-1.) -P-XN/2.-. 25) #1206125 
GO TO 34 
X=-(6.* (R-1.) +P+XN/2.4.25) #120125 
Y=Y+.00014 *X*X 
KCl) = 378 
A(Z) =(Y-.25)/8. 
RETORN 
END 
SUBKOUTINE O'UTLIN (YEAR) 
COMMON A (8) ,22 (9) ,BUP (2048) 
CALL (2L0T (05.50% 53) 
CALLA PLOT (30.5, 14.5, -3) 
GREE PLOT. Oe eo)e 
CALE. PLOT (—24..2,0%, 2) 
CALL PLOT (-5.35,-4.46,3) 
GALE eeLor(-o095, 5-754 2) 
CALL PLOT (-27.91 »-5..25 3) 
CALL PLOT(—21.9'1,6-6,7 2) 
CALE CURCER (—32 9 ly 2OF, 300-705 ¢8050710-0,.020) 
CARLIN PLOT (19 02979-2543) 
GAUL SPLOT (—15.2,. 10.52) 
CALL PLOT (-—13.7,-9.3, 3) 
ChE PEOR (—-13 <7 710.552) 


GALL SYMBOL (-21.,--1,-08, WHHFIFTH MERIDIAN,O., 14) 


ENE SYnBOL ( —23..72 2-13, 10,094,90-,—-1) 
CRUD OBO (G03 07 2, 2.18, « 1eoy 12M G0, 2) 
CALL S¥PRDOL ( —21-20,6.13,-08,003,90.,-1) 
GALL SYP SOL (-2 1.20, 0-18, .08, "CZHROSE*,90.,7) 
CALL SYMBOL (-15.2,6.82,-C8,003,50.,-1) 

CRLL SYMBOL (-15.2,6-27,-08,°STETTLER' /90.,8) 
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END OF FILE 


DS 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALI 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


SYMBOL (-18.2,2.23,-08,003,90.,-1) 

SYMBOL (-18.2,2-28,.08, *PONOKA', 90., 6) 
SYMBOL (-17.9,-1-25, -08,003,90.,-1) 

SYMBOL (-17.9,-1.20, 08, "hIMBEY",90.,6) 
SYMBOL (-15.63,-4-66,.10,004,90.,-1) 
SYMBOL (-15.603,—-4.81,0125,"*RE*,900, 2) 
SYMBOL( -13.91,.61,.10,011,96.,-1) 
SYMBOL (-13-91,.66,-125, "OF*,90.,2) 

SYMBOL (-10.61,-3.34, 068,003, 90.,-1) 
SYMBOL (-10 .61,-3-25, -G8, "SUNDRE*,90.,6) 
SYMLOL (-10.58,-.57, -08,003,90.,-1) 
SYPEOL( -9.84,3.90,.08,003,90.,-1) 

SYMBOL (-9.84,3.95,-08,* THREE FILLS',90.,11) 
SYMBOL( -7.7,7-0,-08,003,50.,-1) 

SYMBOL (-7.7,7.05, 08, "DRUMBZLLER*,90., 10) 
NUMBER (-3-5,11.65,-16,YEAR ,90.,-1) 


CALL SYMBOL (-4.69,-. 12,.10,004,90.,-—1) 


CIMA Sydwsens (NOS -S5O07 pc NZS 4 Ke EEO saa 


RETURN 


END 
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Appendix II 


HAILSWATHS USED IN THE ANALYSIS 
The four hailswaths used in this dissertation are given 
by date in this appendix. The numbers denote the largest 
Size of hail observed by the farmers. The numbers are 


associated with familiar objects and are denoted by: 


Number Hailsize Name Assumed Size Range (cm) 3 
1 Shot ol = 13 

2 Pea o4 - 1.2 

3 Grape Lede =e 2a) 

4 Walnut aot 3.2 

= Golfball 3335 ov ere 

6 Larger than Golfball > Died 


0=is hail of unknown size 
R=rain only 
X=no rain or hail 
*Denotes the hailpad locations 
The vertical line on the figure is the 5th meridian. 
The impact energy densities used in this dissertation for 


each hailswath are: 
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20. July, 1975 
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Figure II.1 The hailifall and hailpad locations for the 24 
June 1974 hailswath 
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.3 The hailfall and hailpad locations for the 3 
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Figure II 


July 1975 hailswath 
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Appendix III 


PERCENTAGE FREQUENCY OF OCCURRENCE OF IMPACT ENERGY DENSITY 


The following pages contain the percentage frequency of 
occurence of impact energy density for the four hailswaths 
examined in this dissertation. The impact energy densities 

2 


are grouped in classes of 25 J m“. The best fitted two- 


parameter gamma distribution is denoted by the "triangles". 
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Figure III.1 The percentage frequency of occurrence of 
impact energy densities for 24 June 1974 hailswath 
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Figure III.3 The percentage frequency of occurrence of 
impact energy densities for 3 July 1975 hailswath 
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Appendix IV 


CUMMULATIVE FREQUENCY OF IMPACT ENERGY DENSITY 


The cumulative frequency of impact energy density for 
the four hailswaths is given in this appendix. The "x" 
represent the original data, while the "triangles" denote 


the best fitted two-parameter gamma distribution. 
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Figure IV.1 The cumulative frequency of impact energy 
density for 24 June 1974 hailswath 
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Figure IV.2 The cumulative frequency of impact energy 
density for 5 July 1974 hailswath 
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Figure IV.3 The cumulative frequency of impact energy 
density for 3 July 1975 hailswath 
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Figure IV.4 The cumulative frequency of impact energy 
density for 20 July 1975 hailswath 


Appendix V 


SAMPLING WITHOUT REPLACEMENT 

For sampling without replacement, the following figures 
are given for the four hailswaths: a) Probability density 
function - here the impact energy density of the sample 
means are grouped in 5 Jm™* classes; b) Probability 
distribution function; c) Standard deviation of the sample 
means as a function of the Sample size; and d) Standard 
deviation of the sample means for the number of hailpads per 


1000 km2, 
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Figure V.1 The probability density function of sampie means 


for 24 June 1974 hailswath 
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Figure V.2 The probability distribution function for 24 
June 1974 hailswath 
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Figure V.3 The standard deviation of the sample means as a 
fandeitdh of the sample size for the 24 June 1974 hailswath 
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Figure V.5 The probability density function of sample means 
for 5 July 1974 hailswath 
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Figure V.6 


200 

JUL OS 74 

WITHOUT REPLACEMENT 
180 N= 1000 


a 
= 
an 


ENERGY DENSITY 


10! 
NUMBER 


aad 
+++ 
x Kx 
eee 
tee 
REX 
Z22 


+t ———+——+__+—_+— 
48 6 7 a8 


154 


the probability, distribution, funchion for 5,.July 


1974 hailswath 


f 


La os) oe parce 4 . dé Me i, Led 
; ' haa a. 
A *= Y ev rh 


155 


o JUL OS 74 
sf WITHGUT REPLACEMENT 
‘i N= 1000 

ros) 


STANDARD DEVIATION 


2 


Figure V.7 The standard deviation of the sample means as a 
function of the sample size for the 5 July 1974 hailswath 
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Figure V.8 The standard deviation of the Sample means for a 
Sample density in 1000 km2 (5 July 1974 hailswath) 
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Figure V.9 The probability density function of Sample means 
for 3 July 1975 hailswath 
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Figure V.11 The standard deviation of the Sample means as a 
function of the sample size for the 3 July 1975 hailswath 
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Figure V.12 The standard deviation of the sample means for a 
sample density in 1000 km2 (3 July 1975 hailswath) 
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Figure V.13..The probability density function of Sample means 
for 20 July 1975 hailswath 
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Figure V.15 The standard deviation of the sample means as a 
function of the sample size (20 July 1975 hailswath) 
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Figure V.16 The standard deviation of the sample means for a 
sample density in 1000 km@ (20 July 1975 hailswath) 
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Appendix VI 


SAMPLING WITH REPLACEMENT 


For sampling with replacement, the following figures 


are given 
function; 
deviation 
size; and 


number of 


for the four hailswaths: a) Probability density 
b) Probability distribution function; c) Standard 
of the sample means as a function of the Sample 
d) Standard deviation of the sample means for the 


hailpads per 1000 kme2, 
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Figure VI.1 The probability density function of sample means 
for 24 June 1974 hailswath 
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Figure VI.2 The probability distribution function for 24 
June 1974 hailswath 
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Figure VI.3 The standard deviation of the sample means aS a 
function of the sample size (24 June 1974 hailswath) 
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Figure VI.4 The standard deviation of the sample means for a 
sample density in 1000 km2 (24 June 1974 hailswath) 
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Figure VI.5 The probability density function of Sample means 
for 5 July 1974 hailswath 
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Figure VI.7 The standard deviation of the sample means as a 
function of the sample size (5 July 1974 hailswath) 
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Figure VI.8 The standard deviation of the sample means for 
a Sample density in 1000 km2 (5 July 1974 hailswath) 
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Figure VI.9 The probability density function of sample 
means for 3 July 1975 hailswath 
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Figure VI.10 The probability distribution function for 3 
July 1975 hailswath 
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Figure VI.11 The standard deviation oté thie Sample means as 
a function of the Sample size (3 July 1975 hailswath) 
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Figure VI.12 The standard deviation of the sample means for 
a Sample density in 1000 km2 (3 July 1975 hailswath) 
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Figure VI.13 The probability density function of sample 
means for 20 July 1975 hailswath 
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Figure VI. 14 
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Figure VI.15 The standard deviation of the sample means as 
a function of the sample size (20 July 1975 hailswath) 
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Figure VI.16 The standard deviation of the sample means for 
a sample density in 1000 km2 (20 July 1975 hailswath) 


Appendix VII 


SAMPLING FROM A FITTED DISTRIBUTION 
For sampling for a fitted two-parameter gamma, the 
following figures are given for the four hailswaths: a) 
Probability density function; b) Probability distribution 
function; c) Standard deviation of the sample means as a 
function of the sample size; and d) Standard deviation of 


the sample means for the number of hailpads per 1000 km2. 
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Figure VII.1 The probability density function of sample 
means for 24 June 1974 hailswath 
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Figure VII.2 
June 1974 hailswath 
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Figure VII.3 The standard deviation of the sample means as 
a function of the sample size (24 June 1974 hailswath) 
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Figure VII.4 The standard deviation of the sample means for 
a sample density in 1000 km? (24 June 1974 hailswath) 
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Figure VII.5 The probability density function of sample 
means for 5 July 1974 hailswath (Impact energy density 
values are grouped in 5 J m~* classes) 
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Figure VII.6 The probability distribution function for 5 
July 1974 hailswath 


ay 


WY 


. f are: st GA <2 f) ri , 


wy a, re Sih 
ia te 
j : ae ed re 


10? 


JUL 05 74 


FITTED DISTRIBUTION 


66. 7-88 


N= 1000 


STANDARD DEVIATION 
10! 


G4 


5 6789 


10° 1Q! 102 
NUMBER 


j lati le means as a 
Figure VII.7 The standard deviation of the sanpl 
function of the sample size (5 July 1974 hailswath) 


4 ines hi Wh 1 
i ae > a) ee 
Miah A, Oy. 7 ae 
i \ ; eae Tr iy, 
mh, ue , "> wii " i" ; 
, voce A | } ij 7 : Am BiG aa 
“ Ph ee 
yay Ae haa: 
y 1 as hy 4@ 
. eo Ue Rene AP yl) 
, Bd OGL nS 8 
oa Le . ed | 
re al & Ht to 


bf 20 JUL. he ae 


WOITUBLATeIO GaTTI3 Bade" A) 
ws ea 
oor == Uae ea 


1960 


10? 


Zs 
oa) ie 
aes 
[— 
‘a 
| ee | 
= 
uJ 
(sa) 
ra Ue 05074 
(ale 
Qe FITTED DISTRIBUTION 
GE 
am 
= N= 1000 
‘ai 
-— 
<p) 
O 


10° 10! 102 
NOWG Era, OOO Kh MaxKM 


Figure VII.8 The standard deviation of the sample means for 
a sample density in 1000 km2 (5 July 1974 hailswath) 
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Figure VII.9 The probability density function of sample 
means for 3 July 1975 hailswath (Impact energy density 
values are grouped in 5 J m~* classes) 
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Figure VII.10 The probability distribution function for 3 
July 1975 hailswath 
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Figure VII.11 The standard deviation of the sample means as 
a function of the sample size (3 July 1975 hailswath) 
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Figure VII.12 The standard deviation of the sample means for 
a sample density in 1000 km? (3 July 1975 hailswath) 
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Figure VII.13 The probability density function of sample 
means for 20 July 1975 hailswath (Impact energy density 
values are grouped in 5 J m™* classes) 
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Figure VII.14 The probability distribution function for 20 
July 1975 hailswath 
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Figure VII.15 The standard deviation of the sample means as 
a function of the sample size (20 July 1975 hailswath) 
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Figure VII.16 The standard deviation of the sample means for 
a sample density in 1000 km2 (20 July 1975 hailswath) 


Appendix VIII 


STRATIFIED SAMPLING WITH REPLACEMENT 


For stratified sampling with replacement, the following 
figures are given for the four hailswaths: a) Probability 
density function; b) Probability distribution function; c) 
Standard deviation of the sample means as a function of the 
sample size; and dj) Standard deviation of the sample means 


for the number of hailpads per 1000 km2. 
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Figure VIII.1 The probability density function of sample 
means for 24 June 1974 hailswath (Impact energy density 
values are grouped in 5 J m~* classes) 
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Figure VIII.2 The probability distribution function for 24 
June 1974 hailiswath 
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Figure VIII.3 The standard deviation of the sample means as 
a function of the sample size (24 June 1974 hailswath) 
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Figure VIII.4 The standard deviation of the sample means for 
a sample density in 1000 km2 (24 June 1974 hailswath) 
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Figure VIII.5 The probability density function of sample 
means for 5 July 1974 hailswath (Impact energy density 
values are grouped in 5 J m~* classes) 
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Figure VIII.6 The probability distribution function for 5 
July 1974 hailswath 
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i VIII.7 The standard deviation of the sample means as 
Agee ze (5 July 1974 hailswath) 
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Figure VIII.8 The standard deviation of the sample means for 
a sample density in 1000 km? (5 July 1974 hailswath) 
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Figure VIII.9 The probability density function of sample 
means for 3 July 1975 hailswath (Impact energy density 
values are grouped in 5 J m~* classes) 
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Figure Vili. 10 *The probability distribution functionsfor 3 
July 1975 hailswath 
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Figure VIII.11 The standard deviation of the sample means as 
a function of the sample size (3 July 1975 hailswath) 
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2 The standard deviation of the sample means 


j VFI of 
Figure (3 July 1975 hailswath) 


for a sample density in 1000 km? 
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Figure VIII.13 The probability density function of sample 
means for 20 July 1975 hailswath (Impact energy density 
values are grouped in 5 J m~? classes) 
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Figure VIII.14 The probability distribution Punct ion» for® 20 
July 1975 hailswath 
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Figure VIII.15 The standard deviation of the sample means as 
a function of the sample size (20 July 1975 hailswath) 
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